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ABSTRACT 

The nature of the low- to intermediate-luminosity (Lx ~ 10 32 34 erg s 1 ) source population 
revealed in hard band (2- lOkeV) X-ray surveys of the Galactic Plane is poorly understood. 
To overcome such problem we cross-correlated the XMM-Newton 3XMM-DR4 survey with 
the infrared 2MASS and GLIMPSE catalogues. We identified reliable X-ray-infrared asso¬ 
ciations for 690 sources. We selected 173 sources having hard X-ray spectra, typical of hard 
X-ray high-mass stars (kT> 5keV), and 517 sources having soft X-ray spectra, typical of 
active coronae. About 18 % of the soft sources are classified in the literature: ~ 91% as stars, 
with a minor fraction of WR stars. Roughly 15 % of the hard sources are classified in the lit¬ 
erature: ~ 68% as high-mass X-ray stars single or in binary systems (WR, Be and HMXBs), 
with a small fraction of G and B stars. We carried out infrared spectroscopic pilot observations 
at the William Herschel Telescope for five hard X-ray sources. Three of them are high-mass 
stars with spectral types WN7-8h, Ofpe/WN9 and Be, and Lx ~ 10 32 -10 33 erg s 1 . One source 
is a colliding-wind binary, while another source is a colliding-wind binary or a supergiant fast 
X-ray transient in quiescence. The Be star is a likely y-Cas system. The nature of the other 
two X-ray sources is uncertain. The distribution of hard X-ray sources in the parameter space 
made of X-ray hardness ratio, infrared colours and X-ray-to-infrared flux ratio suggests that 
many of the unidentified sources are new y-Cas analogues, WRs and low Lx HMXBs. How¬ 
ever, the nature of the X-ray population with K s > 11 and average X-ray-to-infrared flux ratio 
remains unconstrained. 

Key words: Stars: emission-line. Be; Wolf-Rayet - X-rays: stars, binaries. 


1 INTRODUCTION 

The landscape of the Galaxy is dominated by relatively nearby 
stellar coronae at low X-ray luminosities (L x < 10 31 erg_s_*) 

and soft energies (< 2 keV ) (e.g. iMotch et al] 1 19971 : 
iNebot Gomez-Moran et al.l 120131) . In contrast, at high X-ray 
luminosities (L x > 10 35 erg s _1 ), the make-up of the Milky Way is 
dominated by high-mass X-ray binaries (HMXB) and low- mass 
X-ray binaries (LMXB) jGrimm et aD l2002l : lGlifanovll2004l) . On 
the other hand, the nature of the low- to intermediate-luminosity 
source population revealed in hard band (2-lOkeV) X-ray 
surveys of the Galactic Plane, such as those carried out by 
ASCA (fS ugizaki et alj|200lh and more rece ntly by Chandra (e .g. 
lEbisawa et alJl2005l) and XMM-Newton (e.g. lHands et alj[2004h is 
only partly understood. Importantly, in this L x range we expect to 
find systems predicted by evolutionary scenarios of low and high 
mass X-ray binaries. For instance, we have not yet detected the 
long-lived wind-accreting low X-ray luminosity stages preceding 


* Based on service observations made with the WHT operated on the island 
of La Palma by the Isaac Newton Group in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofsica de Canarias. 
t E-mail: ada.nebot@astro.unistra.fr 


or following the bright phase during which they become conspicu¬ 
ous. The common envelope spiral-in creation channel for low-mass 
X-ray binaries predicts the existence of pre-LMXBs which could 
radiate as much as ~ 10 32 erg s _1 in hard X-rays through accretion 
of stellar wind onto the neutron star or its magnetosphere (see e.g. 
iTauris & van den Heuvell 120061) . Likewise, evolution theories of 
high mass X-ray binaries foresee that about 10 6 wind accreting 
binaries containing a main sequence star and a neutron star or a 
black hole populate the Galaxy dPIahl et al jl200^') . Moreover, Be 
+ white dwarf systems are expected to be ten times more frequent 
than Be + neutron stars binaries jRaguzov3l200 ill , whereas none is 
known so far in our Galaxy. With mass accretion rates potentially 
similar to cataclysmic variables (CVs), many of these systems 
should emit copious amounts of hard X-rays. A census of the 
contribution of the different populations to the Galactic hard X-ray 
emission is hence needed if we want to test the validity of these 
evolutionary models. 


Active binaries such as RS CVn systems, and CVs are expec¬ 
ted to contribute si gnificantly to the lo w luminosity hard X-ra y 


source population 

Sazonov et al. 2006;; Revnivtsev et al. 

2009). 

Massive objects such as magnetic OB stars (Gagne et al. 

2011), 

colliding-wind binary O or Wolf Rayet stars (Skinner et al 

20 13 ; 


iMauerhan et alJuOlOl) . y-Cas like objects, being so far the best 
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Be + WD candidates jMotch et~all 120071) . and X-ray transient 
binaries in quiescent state have also been identified in this lu- 
mino s ity range. However, d espite significant efforts I Mot ch et alJ 
120101 : Anderson e t alJI/Olil) . our knowledge of the makeup of the 
Galactic low- to medium-luminosity hard X-ray source population 
associated to massive stars is restricted due to the lack of fully 
identified X-ray source samples. First, the space density of each of 
these sub-classes of X-ray sources is basically unconstrained and, 
second, it remains to be seen whether any new types of Galactic X- 
ray-emitting massive stars, such as those predicted by evolutionary 
theories (e. g. Be + white dwarf systems), are present within cur¬ 
rent hard X-ray Galactic surveys. Owing to sensitivity limits and 
localisation accuracies such studies can only be carried out in our 
Galaxy. 

In this paper we investigate the properties of this important 
class of massive X-ray sources of low- to intermediate-luminosity 
by observing candidates selected from the most recent X-ray 
and infrared surveys. Our target selection, based on the cross¬ 
correlation of the 3XMM-DR4 catalogue with the 2MASS and 
GLIMPSE catalogues, is described in Section [2] As part of a pilot 
study we obtained infrared spectroscopic observations of a small 
sample of hard X-ray sources (see Section[3]i, and produced images 
of the target environment (see Section |4j. which helped to classify 
the sources (see Section^- We discuss our results in Section[6]and 
conclude in Section[7| 


Table 1 . Completeness and spurious fraction of X-ray-infrared associations. 


3XMM-2MASS Pj 


50% 

60% 

70% 

80% 

90% 

Completeness 

89% 

85% 

80% 

72% 

58% 

Spurious fr. 

12% 

10% 

8% 

6% 

3% 



3XMM- 

■GLIMPSE P id 



50% 

60% 

70% 

80% 

90% 

Completeness 

68% 

59% 

49% 

35% 

19% 

Spurious fr. 

24% 

19% 

15% 

10% 

5% 


Notes. f Limited to the GLIMPSE footprint. 


band is ~ 2 x 1CT 15 ergcnr 2 s _1 and source positions have a typ¬ 
ical accuracy better than 3" (90% confidence radius). Among other 
parameters, the catalogue contains information on the count rate, 
some flags describing the quality of the observation, the extension 
of the source, and the hardness ratios (HR), X-ray colours defined 
as: 


HRi = 


C i+ , - C, 
C i+ 1 +Ci’ 


( 1 ) 


where C,- stands for the count rate in the energy band i, where 
i = 0.2-0.5, 0.5-1.0, 1.0-2.0, 2.0-4.5, 4.5-12.0 keV. 


2 TARGET SELECTION 

In this section we briefly describe the cross-matched catalogues and 
the cross-matching procedure. We explain our selection criterion 
and list the X-ray sources chosen for further study at the telescope. 

2.1 Cross-matched catalogues 

2.1.1 3XMM-DR4 

ESA’s X-ray Multi-Mirror observatory XMM-Newton was 
launched in December 1999, and started operations at the beginning 
of 2000 jjansen et al.[|200ll ). It has three telescopes, each equipped 
with one X-ray CCD camera, comprising the European Photon 
Imaging Camera (EPIC). Two of the cameras are MOS_CCD de¬ 
tectors (referred to as MOS-1 and MOS-2 cameras, iTumer et alJ 
l200lh and the third is a pn CCD detector (referred to as EPIC-pn 
camera, [striider et alJ feoQl). While the EPIC-pn camera receives 
all the incident light, a grating disperses about 40% of the light 
going to the EPIC-MOS cameras, diverging it to two Reflection 
Grating Spectrometers (RGS. Iden Herder etal .11200ll) . The field of 
view (FOV) of the telescope is about 30' diameter and the EPIC 
cameras are sensitive in the energy range 0.2-12 keV with a spec¬ 
tral resolution E/AE of about 20-50 and a spatial resolution of 
5" (FWHM). The XMM-Newton Survey Science Centre (SSC) 
compiles catalogues of serendipitous sources detected by the EPIC 
cameras on board the XMM-Newton satellite Jwatson et alJl2009h . 
The 3XMM-DR4 catalogue was released on the 15th September 
20lf] and contains 531 261 detections for 372 728 unique X-ray 
sources. Taking into account overlaps the sky coverage is 794 deg 2 , 
about 1.4% of the sky. The limiting X-ray flux in the 0.2-12 keV 


1 http://xmmssc-www.star.le.ac.uk/Catalogue/3XMM-DR4/ 
UserGuide_xmmcat.html 


2.7.2 GLIMPSE 

The Galactic Legacy Infrared Mid-Plane Survey Extraordinaire 
(GLIMPSE iBeniamin et alJ 12003 ) is a Spitzer Space Telescope 
Mission program that aims to map the inner regions of the Galaxy 
in the infrared. GLIMPS E source catalogue ( I, II + 3D) was 
published in June 2009 JSnitzer Sciencel 2009) through VizieR 
(Ochsenbei n et alJ l2000t> . It provides images and magnitudes for 
over 104 million sources in four bands: 3.6, 4.5, 5.8, and 8.0/jm 
(II, 12, 13 and 14). The catalogue covers Galactic longitudes ±65° 
and latitudes ±1° (up to ±3° in the inner regions of the Galaxy). 
Spatial resolution is < 2", and magnitude limits are about J_5.5, 
15.0, 13.0, and 13.0 in each band respectively ichurchwell et al J 
2009). GLIMPSE covers about 1% of the sky with a typical posi¬ 
tion accuracy of 0.3". 

2.1.3 2MASSPSC 

The Two Micron All Sky Survey (2MASS. Iskrutskie et alJl2006 ) 
used two 1.3 m telescopes located in the northern and southern 
hemispheres equipped with near infrared cameras to obtain simul¬ 
taneous photometry in three bands: J (1.25/tm), H (1.65/im) and 
K (2.17/mi), for the entire sky. The 2MASS point source cata¬ 
logue was released in 2003 and contains about 471 million sources. 
The angular resolution is 2" and the survey completeness limits are 
J = 15.8, H = 15.1 and K s = 14.3 mag. The mean position accuracy 
is ~ 80 mas rms dCutri et aill2003l) . 

2.2 Cross-match procedure 

We cross-matched the 3XMM-DR4 catalogue with the 2MASS and 
GLIMPSE catalogues using the method developed by lPineau et alJ 
J201 it) . In brief this method looks for all the catalogue sources 
around the X-ray source position within a distance lower than 
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Figure 1. Completeness (dashed line) and reliability (solid line) of the 
3XMM-2MASS and 3XMM-GLIMPSE associations as a function of the 
threshold in the probability is shown in black and red respectively. 


3.439cr xc , corresponding to a completeness of 99.7%. Here cr X c 
is the combined positional error of the two catalogues that are be¬ 
ing cross-matched added in quadrature. For all the matches found 
within this distance, i. e for all infrared candidates associated to 
a X-ray source, the method computes the likelihood ratio (LR) 
between the probability that the two sources (X-ray and infrared 
sources) are associated (P id ) and thus have the same position, and 
the probability of finding a spurious association at the same angu¬ 
lar distance and with the same magnitude or brighter than the can¬ 
didate. For that purpose, local densities are calculated for sources 
brighter than the candidate. The LR depends on the probability of 
an X-ray source to have a counterpart in the considered catalogue, 
which depends on the object type considered, its properties and the 
relative fractions of the different populations. To estimate P id , it is 
common to build LR histograms and determine the contribution of 
spurious associations through time-consuming Monte Carlo simu¬ 
lations where positions and errors are randomised. In the method 
developed by IPineau et ali i201 li) . the number of spurious asso¬ 
ciations in each LR bin is estimated, assuming a mean local sur¬ 
face density of X-ray sources, as the sum over all possible associ¬ 
ations of the ratio between the surface element (convolution ellipse 
defined by the positions and their errors) and the field of view area. 
Once the rate of spurious associations is calculated P id can be es¬ 
timated for each LR bin. 

For sources above a given threshold in the identification prob¬ 
ability the method provides us with the completeness and the frac¬ 
tion of possible spurious associations in the cross-matched survey 
(see Fig. |T}. In Table Q] we list the survey completeness and the 
fraction of spurious catalogue identifications as a function of P id . 

In practice, we excluded sources uncorrected for satel¬ 
lite attitude errors (task EPOSCOR) and with positional errors 
(SC_POSERR) in the 3XMM-DR4 catalogue larger than 5". For 
the GLIMPSE catalogue we took into account all sources from 
the v2.0 Archive (v2.0_GLMIIA), i.e. sources detected twice in at 
least one band with S/N > 5 and we limited our sample to sources 
brighter than 14 magnitudes in the 4.5/rm band, corresponding 
to a catalogue completeness of 90% dKobulnickv et alj|2013l) . For 
the 2MASS-PSC catalogue we only considered sources with the 
best photometric quality, Qflag = "AAA’, which corresponds to a 
S/N > 10 and brighter than 15.8 magnitudes in the I band, level at 


which the catalogue is 99% complettQ. We only cross-matched the 
3XMM-DR4 and the 2MASS catalogues in areas also covered by 
GLIMPSE. 

To speed up the process we made use of HEALPix Multi Or¬ 
der Coverage maps (MOC0 which provided us with the coverage 
of a survey within a given minimum pixel resolution. Before per¬ 
forming the cross-match we pre-selected XMM-Newton sources 
within the footprint of the GLIMPSE survey. We accessed cor¬ 
responding MO Cs fr om the interactive software sky atlas Aladin 
(Bonn arel et alj|20 00), and selected sources wi thin the MOC mak¬ 
ing use of the function inMoc in TOPCAlQ jTavlorlf20 1 lh . The 
overlap between 3XMM-DR4 and GLIMPSE footprints covers 
0.1% of the sky, and contains 31 115 unique X-ray sources. 

We found 15 018 possible 3XMM-2MASS pairs (for 9 437 
unique 3XMM sources) and 56 751 3XMM-GLIMPSE possible 
pairs (for 16 668 unique 3XMM sources). In the left panel of Fig. [2] 
we show magnitude histograms of all the possible pairs above a 
certain associated probability. The higher the cut in P, d the brighter 
the sub-sample is, implying that we can only find reliable X-ray- 
infrared associations, i. e. a small fraction of spurious associations, 
for relatively bright infrared sources. 


2.2.1 3XMM, 2MASS, GLIMPSE associations 

Once all the possible 3XMM-2MASS and 3XMM-GLIMPSE asso¬ 
ciations were obtained, the next step was to identify each 3XMM- 
2MASS-GLIMPSE unique association. For that purpose, we pro¬ 
ceeded as follows. 

We limited the sample to X-ray good quality point-like 
sources above a 4a detection. For assuring a 4<r detection we se¬ 
lected sources with parameter sc_det_ml > 10. We defined sources 
as good X-ray detections if sum_flag < 2, where sumjlag contains 
the summary information on 12 flags set automatically and manu¬ 
ally which indicate the quality of the X-ray detection. We restric¬ 
ted our sample to point-like X-ray sources, for which we imposed 
SCJ3XTENT = 0, i.e. no additional source extent model was con¬ 
volved with the point spread function fit to the source count distri¬ 
bution. 

We restricted our study to associations with 3XMM-2MASS 
P id > 0.6, value at which the expected fraction of spurious as¬ 
sociations is about 10% (see Table QJ. For each of these unique 
3XMM-2MASS associations, we looked for all possible 3XMM- 
GLIMPSE counterparts, based on their unique 3XMM-IAUNAME 
identifier, and we considered the right match the one with the 
highest 3XMM-GLIMPSE P id . In this way, we found 2988 pos¬ 
sible 3XMM-2MASS-GLIMPSE associations. Being aware that 
this could eventually lead to wrong associations we created find¬ 
ing charts and visually inspected them to ensure that the chosen 
GLIMPSE counterpart is the correct one. In the right panel of Fig.|2] 
we plot the distance between 3XMM-2MASS and the 3XMM- 
GLIMPSE pairs. There are a few cases where the distance to the 
2MASS association is clearly larger than that of the GLIMPSE as¬ 
sociation and viceversa. Visual inspection in these cases revealed 
that the 2MASS and GLIMPSE sources are possibly not associ¬ 
ated. We prefer to be cautious about these cases and flag them as 


2 http://www.ipac.caltech.edu/2mass/releases/allsky/doc/ 
sec2_2.html 

3 http://www.ivoa.net/documents/MOC/index.html 

4 http://www.star.bris.ac.uk/-mbt/topcat/ 
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Table 2. Infrared counterparts of X-ray sources. 


Col. 

Name 

Contents 






1 

3XMM 

3XMM IAUNAME 

J174347.4-292309 

J180920.3-201857 

J184541.1-025225 

J185210.0+001205 

J190144.5+045914 

2 

RA 

SC_RA 

265.94780 

272.33492 

281.42126 

283.04190 

285.43566 

3 

DEC 

SC.DEC 

-29.38609 

-20.31585 

-2.87375 

0.20143 

4.98748 

4 

POSERR 

SC_POSERR 

1.21 

0.35 

0.44 

0.81 

0.58 

5 

CTS1 

Count rate in 0.2-0.5 keV 

0.00095 

0.00037 

0.00017 

0.00093 

0.00000 

6 

CTS2 

Count rate in 0.5-1.0 keV 

0.00000 

0.00014 

0.00025 

0.00001 

0.00014 

7 

CTS3 

Count rate in 1.0-2.0 keV 

0.00273 

0.00065 

0.00202 

0.00138 

0.00427 

8 

CTS4 

Count rate in 2.0-4.5 keV 

0.00288 

0.00516 

0.01109 

0.00445 

0.01528 

9 

CTS5 

Count rate in 4.5-12.0 keV 

0.01588 

0.00583 

0.03401 

0.00247 

0.01325 

10 

CTS8 

Count rate in 0.2-12.0 keV 

0.02245 

0.01215 

0.04753 

0.00930 

0.03295 

11 

CTS9 

Count rate in 0.5-4.5 keV 

0.00401 

0.00536 

0.01321 

0.00599 

0.01944 

12 

HR1 

Hardness ratio 1 

-0.690 

-0.572 

-0.695 

-0.833 

1.000 

13 

HR2 

Hardness ratio 2 

0.997 

0.963 

0.947 

0.948 

0.971 

14 

HR3 

Hardness ratio 3 

0.557 

0.831 

0.902 

0.751 

0.581 

15 

HR4 

Hardness ratio 4 

0.653 

0.146 

-0.014 

-0.358 

-0.233 

16 

HR1.ERR 

Error in hardness ratio 1 

0.545 

0.215 

0.283 

0.265 

1.346 

17 

HR2.ERR 

Error in hardness ratio 2 

0.158 

0.058 

0.085 

0.096 

0.047 

18 

HR3.ERR 

Error in hardness ratio 3 

0.163 

0.033 

0.021 

0.062 

0.042 

19 

HR4.ERR 

Error in hardness ratio 4 

0.099 

0.038 

0.034 

0.085 

0.044 

20 

2MASS 

2MASS IAUNAME 

17434750-2923098 

18092038-2018572 

18454111-0252257 

18521006+0012073 19014457+0459147 

21 

d_2MASS 

Distance to 2MASS assoc. 

0.4 

0.18 

0.3 

2.24 

0.290 

22 

PROB.2MASS 

Probability of 2MASS assoc. 

0.989 

0.998 

0.994 

0.670 

0.998 

23 

J 

2MASS J magnitude 

11.145 

11.923 

11.562 

11.333 

11.148 

24 

H 

2MASS H magnitude 

7.897 

9.225 

8.700 

9.487 

9.763 

25 

K s 

2MASS K s magnitude 

6.159 

7.719 

7.014 

7.852 

8.841 

26 

eJ 

2MASS J magnitude error 

0.020 

0.027 

0.024 

0.024 

0.022 

27 

eH 

2MASS H magnitude error 

0.042 

0.023 

0.038 

0.028 

0.021 

28 

eK 

2MASS K s magnitude error 

0.017 

0.023 

0.027 

0.027 

0.021 

29 

GLIMPSE 

GLIMPSE IAUNAME 

G359.4078+00.1051 G010.1568-00.3361 G029.7187-00.0316 G033.1949-00.0701 G038.5439-00.0120 

30 

d.GLIMPSE 

Distance to GLIMPSE assoc. 

0.320 

0.34 

0.12 

2.35 

0.110 

31 

PROB.GLIMPSE Probability of GLIMPSE assoc. 

0.998 

0.999 

0.999 

0.966 

0.995 

32 

11 

GLIMPSE 3.6/vm mag. 

4.891 

6.893 


6.537 

7.873 

33 

12 

GLIMPSE 4.5/tm mag. 

4.798 

6.227 

5.346 

5.462 

7.589 

34 

13 

GLIMPSE 5.8//m mag. 

4.152 

5.516 

4.578 

4.781 

7.290 

35 

14 

GLIMPSE 8.0/im mag. 

4.137 

5.322 

4.284 

4.469 

6.981 

36 

ell 

GLIMPSE 3.6jum mag. error 

0.133 

0.064 


0.128 

0.046 

37 

eI2 

GLIMPSE 4.5/rm mag. error 

0.046 

0.086 

0.067 

0.066 

0.043 

38 

eI3 

GLIMPSE 5.8yum mag. error 

0.015 

0.024 

0.024 

0.030 

0.037 

39 

eI4 

GLIMPSE 8.0jum mag. error 

0.019 

0.024 

0.027 

0.021 

0.025 

40 

fA 

X-ray flux in the 2-12 keV 
(ecf- 9.44 X 1(T 12 ) 

2.12e-13 

1.15e-13 

4.49e-13 

8.78e-14 

3.11 e-13 

41 

fk 

Infrared K s flux 

3.18e-09 

7.58e-10 

1.45e-09 

6.69e-10 

2.69e-10 

42 

log (fx/fk) 

X-ray to infrared flux ratio 

-4.25 

-3.86 

-3.53 

-4.01 

-3.00 

43 

44 

X-Type 

NAME 

S for soft or H for hard 

Simbad name 

H 

H 

H 

H 

H 

45 

Class 

Object class 

Star(?) 

WR* 

WR* 

Star(?) 

Be* 

46 

SpT 

Spectral type 

M2I 

WN7-8h 

Ofpe/WN9 

M2I 

09e-B3e 

47 

References 

Spectral type reference 

This paper 

This paper 

This paper 

This paper 

This paper 


Notes.The meaning of each line is given under the column Contents. We estimated the X-ray flux from the EPIC-pn count rate in the 2-12keV band 
using a certain count rate to energy conversion factor ( ecf ). Among others the ecf depends on the combination of camera and filter used, the X-ray emission 
mechanism, and the Galactic absorption. The ecf was computed assuming a thermal Bremsstrahlung model with kT = 5 keV and an absorption of Nh = 10 22 
atoms/cm 2 {ecf ~ 9.44 x 10 -12 count/s /ergcm _2 s _1 ). The choice of the X-ray model and Nh value was done so as to be representative of hard X-ray 
emitting high-mass stars absorbed by the mean Galactic absorption of the sample, calculated from the mean (H-K s ) and assuming sources have no intrinsic 
absorption and a colour (H-K s ) = 0 independent of their spectral type. A different Nh value yields a different ecf, X-ray fluxes are thus to be taken with caution. 
Class and spectral type come from the literature with exception of the five sources studied in this paper. When available the main Simbad name is given. We 
present here an excerpt of the table which will be fully available in the electronic version. 
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Figure 2. Left panel: number of possible 3XMM-2MASS (top) and 3XMM-GLIMPSE (bottom) associations per magnitude bin for Pid> 0.2,0.4,0.6,0.8, 
and 1 (black, dark grey, grey, medium grey and light grey respectively). Right panel: distance between the 3XMM-GLIMPSE possible associations against the 
3XMM-2MASS associations for pairs with a common 3XMM source. 


dubious for any further study. We also checked whether our asso¬ 
ciations between 2MASS and GLIMPSE sources were in agree¬ 
ment with those given in the GLIMPSE catalogue. In ten cases, we 
had different associations and in another 29 cases no association 
was given in the GLIMPSE catalogue. All these 39 cases corres¬ 
ponded to associations already flagged as being dubious based on 
visual inspection of their finding charts. We calculated the distance 
between 2MASS and GLIMPSE sources and found that most of 
the 3XMM-2MASS-GLIMPSE associations with high P u have a 
distance smaller than 1". 


2.3 Selection of hard X-ray high-mass star candidates 

To distinguish sources that are likely active coronae with soft X- 
ray spectra from high-mass stars with hard X-ray spectra we com¬ 
pared the observed X-ray properties of our samp le with modelled 
data. Following the same approach as in lNebot Gomez-Moran et all 
d2013l) we modelled X-ray spectra for different types of objects: 


• Model 1: a young population (70Myr) of active coronae, for 
which we assumed two temperature thermal emission components, 
one componenent (kT| = 0.5 keV) representing the hot plasma and 
a second hotter component (kT 2 = 0.8 keV) dGiidel et ail 1 9971) . 

• Model 2: RS CVn binaries, for which we took as a model val¬ 
ues of WWDra dDemnsev et al.fl993l ). with (kTi, kT 2 ) = (0.2, 2.1) 
and emisivities EM 2 /EM 1 =0.67. 

• Model 3: a population of hard X-ray massive stars. We chose 
one single temperature plasma kT = 5 keV as a representative value 
of a population of hard X-ray massive s tars containing W R stars 
(see e. g. WR 142, OskinovaetaL 20091: ISokal et alj|2010h and y- 
Cas analogues dLones de Qliveiral2007 ). 


We modelled the emission of these three different populations as 
a function of Galactic absorption and we computed the corres¬ 


ponding hardness^ ratios. For more details we refer the reader to 
iNebot Gomez-Moran et al.l d2013l) . 

Among the 3XMM-2MASS-GLIMPSE associations, a large 
fraction have rather large hardness ratio errors, making it difficult 
to differenciate between soft and hard sources. To avoid such prob¬ 
lem we only considered sources in our sample with errors smal¬ 
ler than 0.2 in HR3 and HR4, a value close to the mean error 
in the 3XMM-DR4 catalogue. We finally restricted the sample to 
sources with errors in the HR2 smaller than 1 and that have been 
detected in the EPIC-pn camera, so as to compute fluxes avoiding 
any possible problems related to the merging of the EPIC-MOS 
and EPIC-pn individual detections. We prefered the EPIC-pn cam¬ 
era due to its higher sensitivity compared to the two EPIC-MOS 
cameras. After applying these cuts, we were left with 690 3XMM- 
2MASS-GLIMPSE unique association^ 

In Figure [3] we show the location of all the 690 sources in 
two hardness ratio versus (H-K s ) diagrams. Since the colour (H- 
K s ) remains within 0 to +0.1 from A to K s pect ral types, in¬ 
dependent on the luminosity class dCovev et al.ll2007l) . variations 
in colour mainly reflect variations in Galactic absorption or in¬ 
trinsic emission by circumstellar matter. Under the assumption 
of no intrinsic colour, i. e. (H-K s ) = 0, we calculated the expec¬ 
ted hardness ratios as a function of (H-K s ), using the formula 
N h = 3.5 x 10 22 xE(H-K s ), for the different kind of objects for 
which we modelled the X-ray spectra. We compared the observed 
values to modelled values for active corona, RS CVn binaries and 
hard X-ray massive stars. Among these three type of objects, only 


5 Associations 3XMMJ174505.4-285117-2MASS 17450537-2851166- 
G000.0088+00.1411. 3XMM J175703.9-295921-2MASS 17570406- 
2959219-G000.3666-02.6797 and 3XMM J174912.4-272537- 
2MASS 17491250-2725357-G001.7001 +00.1051 are dubious, the 
distance between the 2MASS and the GLIMPSE possible counterparts is 
larger than 1". 
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Figure 3. X-ray versus infrared colour-colour diagrams for soft (light grey) and hard sources (dark grey) compared to modelled values for active coronae 
(dotted), RS CVn binaries (dashed) and hard X-ray massive stars (solid). Note that we classified as hard sources those having both hardness ratios HR3 and 
HR4 greater than -0.4. Plotted with filled circles are confirmed stars (2 B and 3 G stars) (light grey), HMXBs (red), Be stars (green), and Wolf-Rayet stars 
(blue). The only AGN from the sample is shown as a filled blue upside down triangle and a Nova as a turquoise filled circle. Unidentified hard X-ray sources 
are shown in dark grey, surrounded by a coloured ring when a candidate classification of the source is available: candidate HMXBs (red), candidate YSO 
(yellow), candidate AGB stars (orange), candidate LMXB (magenta). Targets observed at the telescope are shown with a star symbol. See text for more details. 


hard X-ray massive stars are expected to have HR4 above -0.4 at 
any Galactic absorption (see right panel in Fig. [3}. 

We selected sources with hard spectra, defined by HR3 and 
HR4, indicating thermal temperatures greater than 5keV for Nh 
varying from 1 x 10 20 to 3 x 10 22 cnr 2 . The minimum temperat¬ 
ure limit (HR limit) applied to exclude all normal stellar coronal 
emitters is the following: 

HR3 > -0.4 & HR4 > -0.4 (2) 

Selecting sources with HR3 and HR4 greater than -0.4 implies that 
sources are detected in hard X-rays due to either intrinsic hard X- 
ray emission (if (H-K s )~ 0-0.1) or a combination of intrinsic X-ray 
emission and high extinction (if (H-K s )>0.1). Sources with high 
values of HR3, HR4 and (H-K s ) will probably have a high extinc¬ 
tion. 

A total of 173 X-ray sources fulfill the selection criteria 0 
and we thus consider them hard X-ray sources. The remaining 
sources are considered soft X-ray sources in what follows. In Fig. [3] 
we show the location of soft and hard X-ray sources in the X-ray 
versus infrared coulour-colour diagrams. 

We investigated the nature of these 690 sources. We looked 
for possibles entries in SIMBAD making use of XCatDE@, and we 
complemented our search of known sources making use of VizieR. 
In Table[2]we list all the 690 associations. The content of the table 
is as follows. The XMM-Newton IAUNAME is followed by the 
X-ray positions, errors, count rates and hardness ratios. For each 
XMM-Newton source the best 2MASS and GLIMPSE associations 
are given. 2MASS and GLIMPSE IAUNAMES, magnitudes, the 

6 http://xcatdb.unistra.fr/3xmm/ 


distance between the X-ray and the infrared positions and the prob¬ 
ability of identification are listed. The X-ray (2-12keV) and the K s 
fluxes and the X-ray-to-infrared flux ratio are also given. The main 
SIMBAD name, the type of source and the spectral class are given 
when known. We finally give the spectral type reference used for 
the classification of each source. We present here an excerpt of the 
table which will be fully available in the electronic version. 

We estimated the X-ray flux from the EPIC-pn count rate 
in the 2-12keV band and using a count rate to energy conver¬ 
sion factor (ecf). Among others the ecf depends on the com¬ 
bination of camera and filter used, the X- ray em ission mechan¬ 
ism, and the Galactic absorption (see l Mateos et all 20091) . The ecf 
was computed assuming a thermal Bremsstrahlung model with 
kT = 5keV and an absorption of Nh = 10 22 atoms /cm 2 (ecf ~ 
9.44 x 10 -12 count/s/erg cm^s" 1 ). The choice of the X-ray model 
and N h value was done so as to be representative of hard X-ray 
emitting high-mass stars absorbed by the mean Galactic absorp¬ 
tion of our sample, calculated from the mean (H-K s ) and assuming 
sources have no intrinsic absorption and a colour (H-K s ) = 0 inde¬ 
pendent of their spectral type. A different Nh value and/or a dif¬ 
ferent assumed model yields a different ecf, and thus X-ray fluxes 
presented in Table[2]are to be taken with caution. 

As part of a pilot study we were awarded with one 
night at the William Herschel Telescope (see Section [3}. We 
selected targets visible for at least one hour (DEC >-30°, 
15 h < RA < 20 h). We prioritised infrared bright sources and 
went down the list until the end of the observations. We 
could acquire infrared spectroscopy for the following five 
sources: 3XMM J 174347.4-292309, 3XMM J180920.3-201857, 
3XMM1184541.1-025225, 3XMM J185210.0+001205 and 
3XMM 1190144.5+045914. From now on we will refer to these 
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Table 3. Observational technical details. 


Tel. 

Inst. 

Grism 

Slit 

["] 

Spectral 
range [A] 

Spectral 
res. [A] 

Date 

[dd/mm/yyyy] 

WHT 

LIRIS 

j 

0.75 

11690-13530 

3.6 

03/05/2012 

WHT 

LIRIS 

H 

0.75 

15160-17790 

5.2 

03/05/2012 

WHT 

LIRIS 

K 

0.75 

20490-24110 

7.2 

03/05/2012 


sources as Source #1, Source #2, Source #3, Source #4 and 
Source #5 respectively). Finding charts used to identify the 
infrared target at the telescope are shown in Fig. [4] Observations 
and data reduction are described in the next section. 


3 INFRARED SPECTROSCOPY 

WHT/LIRIS. Infrared spectroscopic observations were carried out 
at the William Herschel Telescope (WHT) equipped with the near- 
IR imager/spectrograph LIRIS (Long-slit Intermediate Resolution 
Infrared Spectrograph) on the 3rd of May 2012 (UT) for five of 
our sources. We obtained spectra with gratings J, H and K with a 
slit width of 0.75" (see Tabled- We positioned our targets at two 
nodding points around which we applied small offsets in order to 
subtract the sky and telescope/mirror emissions. The total expos¬ 
ure times were adjusted to obtain a signal to noise of about 100. 
We used standard MIDAS procedures to flat-field correct and ex- 
tract our spec tra, using the optimal extraction method described by 
iHomel dl986h - We took Xenon-Argon arc lamp exposures to cal¬ 
ibrate in wavelength, and obtained spectra of the standard stars 
HIP89218 and HIP93118 at the beginning, middle and end of the 
observations to clean our spectra from telluric lines. Flux calib¬ 
ration was not possible due to a technical problem of the CCD. 
Spectra are shown in Fig. [5] 


We identified spectral lines in the infrared spectra using 


a diversity of references (Morris et al. 1996; Figeretal. 

1997; 

Hanson et al. 1998|; 

Forster Schreiber 2000; Clark & Steele 

2000; 

Steele & Clark 2001 

) and with help of the atomic databases Nisi!/ 


and http://www.pa.uky.edu/~peter/atomic/. 


4 MAGPIS MID-IR AND RADIO IMAGES 

We looked for possible mid-IR and radio emission of our sources 
in the Multi-Array Galactic Plane Imaging Survey (MAGPIS, 
iHelfand et alj|200q) which imaged portions of the first quadrant of 
the Galaxy at 20 cm and 90 cm with the Very Large Array (VLA). 
We retrieved images from the MAGPIS web sitqj, which also in- 
cludes 6 cm VL A, 21 fim Midcourse Space E xperiment ( MSX, 
lEgan et al.l2003l) . and 24 /rm Spitzer MIPSGAL JC arev et al .120091) 
images. In this way we created composite colour-coded images of 
10 arc-minutes around the five targets studied in this paper (see 
Figs. [6TT0t . 


7 http://www.nist.gov/pml/data/asd.cfm 

8 http://third.ucllnl.org/gps/ 


5 SOURCE CLASSIFICATION 

In this section we describe the spectral type and stellar parameters 
for each of the sources observed at the WHT. We estimated the X- 
ray flux in the 0.2-12keV band assuming a thermal model with 
kT = 5keV (see Section l2~Tb and the Galactic Nh we obtained for 
each object. All parameters are given in Tableland some important 
lines with their equivalent width are listed in Table[5] 


5.1 3XMM J174347.4-292309 (Source #1) 

The infrared spectra of Source #1 reveal absorption of many atomic 
lines typical of late K and early M stars, such as Si I, Mg I, Fel, 
KI, All, Nal, Cal and the CO band (see Fig. We used the 
equivalent width (EW) of the CO band-head at 16 187 A and the 
Fe absorption lines at 22 263 A and 22 387 A as indicators of the 
temperature and luminosity class of the star. We compared the 
depth of the CO band head (EW C016187 ^ = 6 ± 1 A) and the 
sum of the EWs of the Fe I 22 263 A and the Fe I at 22 387 A lines 
(EWFe 22263 AA+EWp e 2 2387 A = 4 ± 1 A) with values observed in 
dwarf, giant and supergiant stars iForster Schreibetj|2000i) . These 
values are consistent with a red giant or a red supergiant with effect¬ 
ive temperature close to 3500 K, corresponding to a spectral type 
M2III-I. The EW of the CO-bandhead at 2.3 g m (EW rn „„ = 
41 ± 5 A) is compatible with a supergiant star dNegueruela et alj 
l 20 lch alt hough such high valu es have also been found in some gi¬ 
ant stars IMessineo et alj|2014D . 

In order to construct the SED we complemented the 2MASS 
and GLIMPSE photometric values with information gathered from 
VizieR. A t a distance of 1.6 " from Source #1 there is a WISE 
detection dWright e~ afl l201Ch . WISE J174347.51 -292309.5 ( ICutril 
2012). WISE provides photometry for the whole sky in four band- 
passes centered at 3.4, 4.6, 12 and 22gm (Wl. W2, W3, and W4 
respectively). With magnitudes equal to 5.458 ± 0.038, 4.118 ± 
0.042, 3.982 ± 0.038, 2.927 ± 0.075 in filters Wl, W2. W3, and 
W4 respectively, WISE J174347.51 -292309.5 Wl and W2 mag¬ 
nitudes are not consistent with GLIMPSE II and 12 magnitudes. 
This discrepancy is likely due to the higher spatial resolution of 
the GLIMPSE survey compared to WISE, rather than due to photo¬ 
metric variability. Nearby resolved GLIMPSE sources may indeed 
be blended and contribute to WISE magnitudes. We fitted the SED 
of So urce #1 with two models of spectral types M0 III and M21 
from lCastelli & Kuruczl d2004l) . with (T cff ,log(g)) = (3800 K, 1.34) 
and (T e ff,log(g)) = (3450 K, -0.06) respectively. The (H-K s ) colour 
indicates a large extinction, 22<Av<27. We fixed extinction to 
A v = 26 magnitudes, equal to the total interstellar extinction in the 
area dSchultheis et al J1 1 9991). According to the empirical spectral 
type radius relation from lvan Belle et all dl999l) . assuming a giant 
star, the radius is 40/ 7 Ro, for an error in the spectral type determ¬ 
ination of ± 1 sub-type. From the fitting flux scaling factor we es¬ 
timated a distance to the star of 310/™ pc. We are not aware of 
the existence of such a nearby high extinction region and we thus 
conclude that the star must be a supergia nt. Supergiant stars ex- 
hibit a wide range of luminosities (see e.g. lMevnet & Maederl200(i 
iMarigo et alj|2008f) and as a res ult a wide range of ab solute mag¬ 
nitudes is expected and observed dLevesaue et alj2 005). Therefore, 
no reliable distance may be derived from the photometry. The high 
extinction observe d indicates tha t the source must be rather dis¬ 
tant. According to iMarshall et aid d2006h the Ak extinction in the 
line of sight (///,/>//) = (359.5° ,0°) reaches a value close to 2.5 at 
a distance larger than 4.5 kpc. At a distance of 4.5 kpc the source 
absolute magnitude Mk is equal to -10.07, a value consistent with 
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3XMM J 174347.4-292309 

» 



3XMMJ 180^03-201857 4fj 


3XMMJ184541.1-025225 



3XMM J185210.0+00120 £ 


© 


Figure 4. Two arc-minute square 2MASS J-band image centred around the five observed targets at the WHT. Positions and 3cr errors are shown with circles, 
in white for 3XMM, in green for 2MASS and in red for GLIMPSE detections. Following the standard convention. North is up and East is left. 



3XMM J! 84541.1 -025225 







WAwl j 







3XMM J190144.5+045914 




t* v w w 




Figure 5. Observed SED (coloured dots) compared to the best fitting absorbed model or Rayleigh Jeans distribution (solid or dashed grey line), followed by 
the WHT infrared spectra in the J, H and K bands. From top to bottom: Source #1, Source #2, Source #3, Source #4, and Source #5. 
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Table 4. Parameters 


Name 

3XMMJ 

SpT 

T c ff 

[K] 

Rad 

[Rel 

A V 

N h 

[xlO 22 cm~ 2 ] 

H 

K s 

M k 

d 

[kpc] 

Count rate 
[cts s -1 ] 

L x (0.2-12 keV) 
[10 32 ergs~'] 

Nature 

174347.4-292309 

M21 

3500 

640 

~26 

4.6 

7.897 

6.159 

-10.2* 

>4.5 

0.022 ± 0.006 

> 6 

? 

180920.3-201857 

WN8 

40 000 

30 ±10 

23 ±3 

4.1 ± 0.5 

9.225 

7.719 

-6.79+ —9.29 

1.6-5 

0.012 ± 0.006 

0.4-4 

WCB 

184541.1-025225 

Ofpe/WN9 

20 000 

79 

26 ±3 

4.6 ± 0.5 

8.700 

7.014 

-7.9 ±0.1 

2.4 ±0.3 

0.047 ± 0.007 

2.3-5.3 

WCB/SFXT 

185210.0+001205 

M21 

3 500 

430 

~ 15 

5.4 ± 0.2 

9.487 

7.852 

-9.06 

11 

0.009 ± 0.002 

< 15 

? 

190144.5+045914 

Be 

- 

- 

15 ±2 

2.7 ± 0.2 

9.763 

8.841 

-3.15--4.43 

1.0-1.9 

0.033 ± 0.002 

0.9-3.0 

-y-Cas 


Notes. We list the 3XMMIAUNAME of the five targets observed at the telescope, the spectral type of the infrared association derived from the spectra taken at 
the WHT, eff ective temperature a nd radius derived from the SED fit, Ay estimated from the infrared colour (see text for more details), Nh calculated using the 
relation from lPred ehl & Schmitt) 119951) . Obse rved H and K s 2MA SS magnitudes are followed by the absolute magnitude in the K s band. We adopted absolute 
magnitudes from [Crowther et a i. IlSIt and lCovev et al.l J2007l) +. Distances are calculated from the adopted absolute magnitudes. Finaly we list the count 
rate and X-ray luminosity in the 0.2-12 keV band, followed by the likely nature of the X-ray source. 


those of red supergiant stars jLevesaue et al.ll2005h . Fixing the dis¬ 
tance, from the best SED fit, we obtained a radius of about 640 R 0 , 
a value consistent with known red supergiant stars. A larger dis¬ 
tance to the source would imply a brighter absolute magnitude and 
a larger radius. 

Source #1 was detected in two different observations with 
a mean count rate of 0.022 ± 0.006 ctss -1 in the 0.2-12keV 
band. We derived an unabsorbed X-ray luminosity of Lx > 6 x 
10 32 erg s~*(N H =4.6 x 10 22 cm' 2 ). No intrinsic absorption was as¬ 
sumed. 

On five occasions, Chandra detected the hard source 
CXOGC 174347.4-292309 (RA = 1 7 :43:47 .5 0, DEC = - 
29:23:09.9, r95%=0.7") d \11111 o et a 0 120091 : Mauerhan et alj 
120091 : iHong et all l2009h at a distance of 1.7" from the XMM- 
Newton position. Given the positional errors of Source #1 (1.2") 
and of CXOGC 174347.4-292309, the Chandra and XMM detec¬ 
tions are likely to be associated with the same source. Importantly, 
the improved Chandra positio n is consistent with the GLIMPSE 
candidate l lMuno et alj l2009h . Long-term X-ray flu x variations 
were detected between different Chandra observations 1 M uno et alj 
i2009}l . The ratio between maximum and minimum observed fluxes 
is in the range 2.3 to 8.2, and has an associated time span of ap¬ 
proximately five years. Maximum flux was detected in September 
2006. The two XMM-Newton observations in which the source 
was detected span less than one year, both of them in 2006 (Feb¬ 
ruary and September). No significant flux variation was detected 
between these two observations nor within the observations. 
Assuming the X-ray luminosity of L x = 6 x 10 32 ergs~* (derived 
from the XMM-Newton observations) as a maximum value the 
ratio between the minimum to maximum flux observed by Chandra 
implies a minimum X-ray luminosity of L x ~ 7.3 x 10 31 erg s _1 . 

5.2 3XMM J180920.3-201857 (Source #2) 

The WHT J, H and K band spectra of Source #2 reveal very broad 
emission lines, in particular from the complex at 12 788 - 12 818 A, 
formed by C III, He I, He II, and HI (see Fig. 0, and from the 
complex at 21 652-21 663 A, formed by Hell, He I and HI. The 
prese nce of such emission lines is typical of WN stars dShara et~aT] 
12009ll . The strong emission He IIT21 891 A line indicates a WN7- 
8h spectral type. We estimated the ratio of the EW of the 21 891 A 
Hell line to the EW of the complex at 21 663 A to be 0.1, and 
between the 21 891A and the 21 126 A Hell lines to be 0.5. These 



Figure 6. Ten arc-minute square RGB image centred around Source #1 (red 
= 24gm, green = 8 fjm GLIMPSE, blue = EPIC-pn XMM-Newton). North 
is up and East is left. 


Table 5. Equivalent widths (in A) of some important emission lines (in /./m). 


Name 

3XMMJ 

HeII+P/3 HeH 
1.281 1.692 

He I He I He I/NIII He II+Bry He H 
1.700 2.058 2.115 2.165 2.189 

180920.3-201857 

71 

2 

9 

19 

64 

8 

184541.1-025225 

90 

- 

3 

17 

40 

2 

190144.5+045914 

6 

- 

1 

9 

10 

3 


values are in agreement with a WN8 spectral type dFiger et alj 

ll997l : lRosslowe & Crowthei)l2015t). _ 

We followed the method in lCrowther et alj J2OO6) to determ¬ 
ine the distance to the source. From the 2MASS c olours and as- 
suming intrinsic colours of WN7-8h stars from ICrowther et alj 
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Figure 7. Ten arc-minute square RGB image centred around Source #2 
(red = radio 20-cm, green = 13-GLIMPSE filter, blue = EPIC-pn XMM- 
Newton). North is up and East is left. 



Figure 8. Ten arc-minute square RGB image centred around Source #3 (red 
= radio 20 cm, green = 24pm, blue = EPIC-pn XMM-Newton). North is up 
and East is left. 


( 2006 1 we calculated the colour excess E(J-K S ) and E(H-K S ). 
We used the extinction relations A K =0.6 7 + Siii lx E(J -KJ and 
Ak = 1.82!““xE(H-K s ) from [indebetouw et al. (20051) to derive 
a mean value Ak= 2.63lQ2g. Using the extinction relation Ak = 
0.114xA v for R v = 3.1 from Cardellietah d 19891) we estimated 
A v = 23 ± 3. iRosslowe & Crowtherl d2015 ) revised the infrared 
absolute magnitudes of the largest sample of WR stars to date. 
For WN7ha and WN9ha stars they obtained a mean value of Mk 
of -7.24 and -6.34 respectively, with a spread of about 0.3 mag¬ 
nitudes at each subtype. From the observed K s magnitude, the mean 
Ak and assuming an absolute magnitude Mk equal to the mean 
value for WN7a and WN9ha stars (M K = -6.79 ± 0.42) we es¬ 
timated a distance to the source of 2.4^?S kpc. 

We compared the observed infr ared spectra with model spec¬ 
tra from lHamann & Grafener ( 2004 ) and found a best match for a 
WNL with an effective temperature of about 40 000 K. Using this 
T cff and a distance to 2.4 kpc implies a radius of 30 ± 10 Rq. Fix¬ 
ing the extinction and distance, no black body can fit the observed 
SED, even when discarding the GLIMPSE II and 12 magnitudes 
which could be confused with surrounding nebular emission. We 
therefore left the extinction as a free parameter and obtained a best 
fit for Ak ~ 3.5 (Ay ~ 30). Assuming a mean absolute magnitude 
for WN7a and WN9ha we obtained a distance « 1.6 kpc, which is 
an extremely small distance for such a high extinction. 

In Figure[7]we show a 10 arc minute composite colour image 
around Source #2. The observed infrared nebulosity is likely as¬ 
sociated with Source #2, i.e. the strong wind and ionising photons 
produced by this hot massive star blow away the surrounding ma¬ 
terial forming a cavity around the star, and ionising the gas in the 
nebulae, which is also seen as radio emission (see Figure[7). Indeed, 
the X-ray source isjn the infrared dark cloud SDCG10.156-0.340 
(IRDC, |Peretto^Fuller 20091). e mbedded in an infrared bubble 
(N 1 from IChurchwell et alj2006l : lDeharveng et alj2010h . Accord¬ 
ing to lDeharveng et all d2012l) this bubble is in an HII region that 
forms part of a more complex structure, the star-forming region 


W31 at spectrophotometric distance 3.4 ± 0.3 kpc d B1 um et al] 
hoQli) . At a distance of about 1.5' from Source #2 the source 
3XMMJ180926.9-201930 was detected by the XMM-Newton as 
an extended X-ray source. This extended source has a very hard 
spectrum and is likely associated to the OB cluster W31 IR Cluster. 
[Stead & Hoarei <[201 I T estimated the cluster to be formed by about 
35 candidates with extinction Aj = 4.3*Jg (A K = l-7 0, ® 2 ). 

Although the lower limit of the extinction Ak we obtained for 
Source #2 is consistent with the upper limit of such a value, the 
distance to the cluster is higher than to Source #2 by at least 500 
pc. Moreover, in a recent study ISanna et al.1 1 20141) find a distance 
to the cluster of 5.0 ± 0.5 kpc based on parallax measurements, a 
value inconsistent with our distance estimate for Source #2. It is 
a bit surprising that Source #2 has a higher extinction than W31 
while being foreground to W31. Assuming that the source is as¬ 
sociated with the cluster, i.e. located at 5 kpc and is absorbed by 
A k ~ 1.7 implies M K ~ -8.39. Under these conditions, no black 
body can fit the observed SED. Leaving the extinction as a free 
parameter yields Ak ~ 3.5 and Mk ~ -9.29, a value higher than 
observed in typical WRs. 

Source #2 was detected in 13 different observations with a 
mean count rate of 0.012 ± 0.006 ctss -1 and with no signific¬ 
ant flux variation. Ass uming A K = 2.63; f ° 33 , c orresponding to 
N h ~ 5.5 x 10 22 cnr 2 dPredehl & Schmittlll995t) . the unabsorbed 
X-ray luminosity is in the range 4 x 10 31 - 4 x 10 32 erg s -1 depend¬ 
ing on the assumed distance (1.6 kpc-5 kpc) to the source. The 
bolometric lumi nosity log(L|,ni /L w ) is in the range of 6.29 to 7.29 
(BCk = -4.2 dMauerhan et alJuOlCh ) yielding log(Lx/Lb 0 i) = 
-8.2, a value lower than the mean value for WN stars but within 
the observed range (Mauerhan et al. l20ld) . 


5.3 3XMM J184541.1-025225 (Source #3) 

The infrared spectra of Source #3 reveal strong and narrow emis¬ 
sion lines of H I, blended with He I, and of some ionised species, in 
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particular Fe II and Mg II. The H band spectra displayed the forbid¬ 
den line [Fell] at 16770 A blended with the HI (11 -4) line. The 
absence of the 16 923 A He II line implies a spectral type later than 
O. The spectra also shows Na 122 060 A and 22 090 A in emission. 
These lines are typ i cal of Ofpe/WN9 a nd early Be and B[e] stars 
dMonis et alJll996l : lHanson et alll 1998h . Whether the source is a 
luminous blue variable star (LBV), i.e. a source displaying spectro¬ 
scopic variations between both spectral types, is an open question, 
and more observations are required to elucidate the nature of the 
source. 

The 90% error circle of Source #3 also contains 
WISEJ18454111-025225.5 which is probably associated with 
the GLIMPSE identification. In Figure [5] we present the observed 
SED. We determined the mean extinction following the same 
procedure as for Source #2. The mean Ak value is 3.0^ 3 . We as¬ 
sumed an effective temperature of 20 000 K and a radius of 79 R e , 
corresponding to the Ofpe/WN9 star CXOGCJ174516.1-284909 
dMauerhan et alj2010n and fitted a black body curve. From the best 
fit we obtained a distance to the source of 2.4 ± 0.3 kpc, implying 
an absolute magnitude M K = -7.9 ± 0. l. lNegueruela et al.l d201 II ) 
studied in detail several nearby sight-lines. These areas showed 
moderate extinction up to distances of 3 to 3.5 kpc and with an 
abrupt increase between 6 to 7 kpc (A K > 1.5) due to dark clouds 
associated to the Scutum arm. At the calculated Ak the source 
would be even more distant. Assuming the source is at about 
2.4 kpc it would be in the inter-arm region. Nevertheless, it could 
be that the source has some local obscuration and or emission, 
especially since the source has [W1]-[W4] = 2.6, a colour atypical 
of interstellar absorption. 

About 3.6" from Source #3 there is the radio emitter 
GPS5 029.718-0.031, detected in the in the 5 GHz VLA survey of 
the galactic plane l iBecker et al.ll 1994h . The radio position and flux 
of GPS5 029.718-0.031 is consistent w ith those of the c ompact 
radio source [WBH2005] 29.719-0.031 dWhite et alj|2005h . about 
l"away from the X-ray position but within the combined position 
errors. The GLIMPSE sourc e is probably asso ciated to the radio 
source GPS5 029.718-0.031 dHoare et alll2012h . The radio emis¬ 
sion might be ascribed to a massive young stellar object (MYSO) 
associated to collimated jets and winds, or it could originate in a 
new born HII region. The detected radio emission at 6 and 21 cm 
(Becker et all l 19941) is consistent with a flat spectral index a ~ 0 
(S r oc v"), indicating a non-thermal emission, and excluding a pos¬ 
sible MYSO nature of the source, objects with a typical a value 
~ +0.6. X-ray emission associated to high-mass stars has been 
sugge sted as an ionising source of UCHII regions dHofner et alj 
120021) . We conclude that the radio emission detected in the com¬ 
bined X-ray plus radio 90% confidence error circle of Source #3 is 
associated to a radio loud ultra compact HII region (UCH II). This 
ultra-compact HII region, detected at radio wavelengths could ex¬ 
plain the observed [W1]-[W4] colours. We adopted a distance to 
the source of 2.4 kpc. 

The three X-ray detections of Source #3 have a mean count 
rate equal to 0.047 ± 0.007 cts s _1 in the 0.2 - 12 keV band, imply¬ 
ing an unabsorbed X-ray luminosity L x = 3.7 + | ^ x 10 32 ergs~‘ . 
We assumed the bolometric correction from Mau erhan et al.l OOlOl ) 
derived for the Ofpe/WN9 star CXOGC J174516.1-284909 and ob¬ 
tained log L bo [ = 6.21 L 0 , and thus log(L x /L bo i) = -7.2, value 
consistent with the canonical value of -7. 



Figure 9. Ten arc-minute square RGB image centred around Source #4 (red 
= MIPSGAL 24 fim, green = GLIMPSE 8.0/rm, blue = EPIC-pn XMM- 
Newton). North is up and East is left. 


5.4 3XMM J185210.0+001205 (Source #4) 

As for Source #1, the infrared spectra of Source #4 reveal absorp¬ 
tion lines of many atomic species typical of late K and early M stars 
(Si I, Mg I, Fel, KI, All, Nal, Cal and the CO band, see Fig.0. 
Besides the Nal 12683 A emission line in the I band spectrum, no 
other emission line could be identified in the spectra. Following the 
same approach as for Source #1 we measured the EW of several 
absorption lines to determine the spectral type and luminosity class 
of the star. With an EW C0 232 27 a = 20 ± 1 A, an Fel 22263 A to 
Fe I at 22 387 A EW ratio equal to 5 ± 2 and EWco 2 . 3 fim = 43 ± 5 A 
the spectrum of the star is consistent with that of a red super¬ 
giant with effective te mperature close to 3 500 K, corresponding 
to a s pectral type M21 dForster Schrei ber 2000; Nes ueruela et alJ 

l20Ki|) . 

Source #4 has the optical counterpart U SNO 0902-0385469 
in the USNO-B 1.0 catalogue dMonet et alj|2003l) . at a distance 
of 0.26" and mean magnitudes equal to 20.9, 17.7, and 15.2 
in the B. R and I filters respectively. The source was detected 
in the 2XMMi-DR3 and was already listed as having IR and 
optical counterparts by INebot Gomez-Moran et all d2013l) . At a 
distance of about 2" from Source #4 there is the WISE detec¬ 
tion WISE J185210.06+001207.3 with magnitudes 6.411 ± 0.038, 
5.404 ± 0.025, 4.675 ± 0.017, and 3.923 ± 0.040 in filters Wl, 
W2, W3 and W4 respectively. WISE magnitudes are in agreement 
with GLIMPSE magnitudes. We constructed the SED and con¬ 
cluded that it cannot be explained by asingle atmosphere compon¬ 
ent when comparing with the lCastelli & Kuruca |2004j) model with 
T e ff = 3450 K and log(g) = -0.06, corresponding to a spectral type 
M2I. Assuming an intrinsi c col our (I-K s ) of 0.87 equal to the val¬ 
ues found for M21 stars ( iNegueruela et al.l20ld) we derived an ex¬ 
tinction Ak~ 1.75 (Ay ~ 151. According to iMarshall et alJ2006h 
the extinction in the line of sight (lII,b II) = (33.25°,0°) reaches 
that value at about 11 kpc, giving us an upper limit on the distance 
to the source. No single atmosphere model is able to reproduce the 
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Figure 10. Ten arc-minute square RGB image centred around Source #5 
(red = MIPSGAL 24pm, green = GLIMPSE 8.1)pm, blue = EPIC-pn 
XMM-Newton). North is up and East is left. 


observed SED from the B to the W4 bands. To match the blue part 
of the SED (B to K s bands) the needed attenuation is Av~ 18, 
displaying an excess at longer wavelengths. To match the infrared 
part of the SED (H to W4 bands) requires Ay ~ 28 displaying an 
excess at shorter wavelengths. The extinction inferred from the in¬ 
frared colours is closer to the extinction obtained through SED fit¬ 
ting in the blue part of the spe ctra than to the red part. Source #4 
is in a dark cloud |Peretto_&FuUei 200S) within a molecular cloud 
and an HII region (Anderson et al. 2008 ), which is consistent with 
a high extinction. It could be that the GLIMPSE source is confused 
with this HII region at longer wavelengths, explaining the disagree¬ 
ment between the observed SED and that of a red supergiant star. 
Assuming an upper limit of 11 kpc to the source and an extinction 
of Ak ~ 1.75 we obtained an absolute magnitude of Mj = -8.2. 
From the SED best fit to the blue part of the spectra we estimated 
that the radius of the source is around 430 R 0 . 

Source #4 was detected in two XMM-Newton observations, 
with mean count rate 0.009 ± 0.002 cts s' 1 , implying an X-ray lu¬ 
minosity of 1.5 x 10 33 erg s' 1 (at a distance of 11 kpc. No variability 
was observed between the two observations. 


5.5 3XMM J190144.5+045914 (Source #5) 

The infrared spectra of Source #5 show the HI Brackett and Pfund 
series in emission, which together with the He I also seen in emis¬ 
sion, indicate that the sourc e is probably an early type Be star, 
with spectral type 09e-B3e t jciark & Stee 1 el 2 00(j| : l Steele & Clarkl 
l200ll) . where the X-ray and UV photons originating from the hot B 
star ionise the circumstellar disc giving rise to strong hydrogen re¬ 
combination emission lines. Photospheric absorption lines are vis¬ 
ible in the K band spectrum. Emission lines of He II, and Fe I are 
also present in the spectra (see Fig. [5}. 

The GLIMPSE and MIPSGAL images at 8 pm and 24/jm 
respectively revealed an extended dark cloud in the direction of 
Source #5 (see RGB composite image in Fig. llOt . At 8 pm, emis¬ 


sion is dominated by poly-cyclic aromatic hydrocarbons, and at 
24pm emission is dominated by very small dust grains. The emis¬ 
sion from this warm dust can explain the high extinction observed 
for Source #5, for which we derived A K = 1 •69^' 3 j (A v = 15 ± 2) 
magnitudes from the (H- K s ) co lour and assuming the intrinsic col¬ 
ours of 09.5 stars from iMartins & Pled J20061 ) (note that these 
colours are independent of the stellar luminosity class). How¬ 
ever, the source was also detected in the A11WISE catalogue, All- 
WISEJ190144.57+045914.8, with magnitudes that are not consist¬ 
ent with a reddened Rayleigh distribution (see Fig. [5}, suggesting 
some local excess, probably from the circumstellar Be disk. There¬ 
fore we cannot use the observed infrared colour as a reliable in¬ 
dicator of the extinction and the Ay value should be considered as 
an upper limit to the extinction. At shorter wavelengths the USNO 
photographic magnitudes are not accurate enough to constrain the 
reddening. 

Late-O and early-B supergiants show absorption lines in the 
infrared spectra, while the bright est blue supergiants show spectra 
with strong wind emission lines (Hanson et al. 120051) . both incon¬ 
sistent with our infrared spectra. Since we don’t know the luminos¬ 
ity class of the star we assumed different tabulated magnitudes for 
09.5 stars with luminosity classes V and III from IMartins & Plea 
( 2006 ) and derived maximum distances of about 1.0 to 1.9 kpc de¬ 
pending on the luminosity class. The possible link of this system to 
the class of y-Cas analogues will be discussed in Section R+Tl 

X-ray emission from Source #5 was fi rst detected by ASCA, 
AX J190144+0459 iSugizaki et alj Eooi ). as extended diffuse 
emission, however, the extension was n ot co nfirmed by subsequent 
XMM-Newton observationj jfamaguchi_etji p2004l) . A power law 
model with T= 1.7 lYamaguchi et a iJ ( 12004} ) derived Nh = 3.3 x 
10 22 cm' 2 and a N H corrected flux of 5.1 x 10' 13 erg ernes' 1 in 
the 0.5-10keV band. In the range of distances we derived this im¬ 
plies L x ~ 6.4 x 10 31 - 2.1 x 10 32 erg s' 1 . ICorral et al] J2014I) 
fitted the source to different models and present the results in the 
XMM-Newton spectral-fit catalogu^ In the energy range from 
0.5 to 10 keV, the best fit is found for an absorbed thermal 


- 2.23 keV, N H = 2.99^28 x 10 2 


atoms cm 


plasma with kT = 9.27 
with a corresponding absorption corrected flux of 3.14 ± 0.34 x 
10' 13 erg ernes' 1 , implying L x in the range of ~ 4.0 x 10 31 to 
1.3 x 10 32 ergs' 1 . 


6 DISCUSSION 

In this section we discuss the possible origin of the hard X-ray 
emission observed in these five X-ray sources, where infrared coun¬ 
terparts pointed to three high-mass stars and two descendants of 
high-mass stars. We then discuss the nature of all the hard X-ray 
sources found in this study. 


6.1 High-mass stars 

The origin of hard X-ray emission in massive stars is still in debate. 
There are mainly three competing scenarios: intrinsic X-ray emis¬ 
sion, accretion on a compact object and colliding winds in a binary 
system. These three scenarios have been discussed in great detail in 
iMauerhan et al.1 l20ld) . We here briefly summarise the main prop¬ 
erties of systems in the different scenarios. 

The intrinsic X-ray emission of early type stars is typically 


9 http://xraygroup.astro.noa.gr/Webpage-prodec/index.html 
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thermal, soft (kT ~ 0.5 keV) and due to_radiative shocks occurring 
within the high velocity wind dLucv & White|[l980l) . This mech¬ 
anism is thought to be at work in_a wide range of massive stars, 
including Wolf-Rayet stars dSkinner et all l2Q10l) . Unusually hard 
X-ray emission, although rare in single objects, has been detected 
in some massive stars (e.g. S'-Orionis C). As a possible explanation 
to the hard X-ray emission (kT ~ 3 keV) it has been suggested that 
the high magnetic fields observed (> 1 kG) could confine and chan¬ 
nel the stellar wind from the poles into the equator where a shock 
could heat the plasma up to 10 7 K, thus producin g hard X-ray emis¬ 
sion dGasne et alj2005l) . lMauerhan et al,l d2010h calculated that the 
strength of the magnetic field should be > 5 kG in order to con¬ 
fine the fast winds observed in the WR and O stars of their sample, 
which is similar to ours. Since so far there is no observational evid¬ 
ence of such high magnetic fields in WR stars, we disregard this 
scenario for our WR candidates. 

In the second scenario, a compact object accretes matter from 
a massive star, i.e. high-mass X-ray binaries (HMXB). There are 
three types of HMXBs: those accreting material from the circums- 
tellar decretion disc of a Be star (Be-HMXB), those accreting ma¬ 
terial from the wind of a supergiant star (SGXB), and those accret¬ 
ing from a Roche lobe filling supergiant star ( IChatvl[201 ill . Since 
Be-HMXBs display L x > 10 33 erg s -1 , the low observed X-ray lu¬ 
minosity of the Be star counterpart of Source #5 indicates that the 
source is most probably not a classical Be-HMXB in quiescence. 

Supergiant high-mass X-ray binaries (SGXBs), i.e. a neutron 
star or a black hole accreting from the strong radiative wind emitted 
by an OB supergiant star or by a WR star, have persistent X-ray 
luminosities L x ~ 10 35 - 10 36 ergs -1 , show rapid varia tions and 
are rather constant on long time scales (see lSeuera et alfeood . and 
references therein). The X-ray luminosities observed in our sample 
of massive stars are too low to be compatible with SGXBs. 

However, a sub-class of accreting sources associated with su¬ 
pergiant stars, the Supergiant fast X-ray transients (SFXTs) show 
quiescent luminosities in the range L x ~ 10 32 - 10 33 ergs -1 , con¬ 
sistent with those of our five stars, while exhibiting flares last- 
ing less than one day and reaching up to L x ~ 10 36 erg s -1 

iNegueruela et all200a) . The physical origin of the energetic flares 
is still not clear, it can be either due to the accretion of clumps 
from the supergiant wind or to some gating mechanism acting 
close to the compact accreting object. The X-ray light curves of 
the only source possibly associated with an early supergiant coun¬ 
terpart (Source #3) does not show any flaring activity characteristic 
of SFXTs. However, the wide range of inactivity duty cycles dis¬ 
played by SFXTs l lRomano et 10 SQM) does not allow us to ex¬ 
clude a SFXT nature for this particular source. 

Likewise, none of the five sources reported in this paper may 
qualify as an HMXBs undergoing Roche lobe overflow which usu¬ 
ally exhibit relatively steady X-ray luminosities of the order of 
L x ~ io 36-38 erg s -1 . 

The third possible explanation for hard X-ray emission in 
high-mass stars is the colliding-wind binary interpretation. In this 
scenario opposing winds generated in two massive stars collide 
producing a hot shock that is the source of hard thermal X-rays. 
There are many cases reported in the literature of such systems 
iMauerhan et all 20 id) . where WR + O stars are the most frequent 
association. X-ray luminosities are in the 10 32 - 1(L 4 erg s -1 range, 
i.e. compatible with the X-ray luminosities observed in Source #2 
and Source #3. Although in our spectra we did not detect any addi¬ 
tional spectral lines hinting at the presence of a O companion star, 
the lower optical luminosity of these objects would not allow us 
to detect the OB dwarf spectrum. Source #2 has very broad emis- 
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Figure 11. Infrared K s magnitude versus X-ray-to-infrared flux ratio (in the 
2-12 keV band) for hard X-ray sources. Plotted with filled circles are con¬ 
firmed stars (2 B and 3 G stars) (light grey), HMXBs (red). Be stars (green), 
and Wolf-Rayet stars (blue). The only AGN from the sample is shown as 
a blue filled upside down triangle and a Nova as a turquoise filled circle. 
Unidentified sources are shown in dark grey, surrounded by a coloured ring 
when a candidate classification of the source is available: candidate HMXBs 
(red), candidate YSO (yellow), candidate AGB stars (orange), candidate 
LMXB (magenta). Targets observed at the telescope are highlighted with 
a star symbol colour coded according to their classification. The effect of 
extinction assuming a population of hard X-ray massive stars (model 3 in 
Scction [2~3l is shown with an arrow. 


sion lines, associated to a very fast wind. It could also be that these 
emission lines are blended with a second WR star, which would not 
be resolved given the available resolution. 

Finally, y-Cas analog ues con stitute a distinct class of hard X- 
ray emitting Be stars IlSmith et alJl2015n . Confined into a narrow 
range of spectral types ranging from BOe to B1.5e with dwarf or 
giant luminosity classes, their X-ray emission appears essentially 
thermal with temperatures larger than lOkeV and luminosities of 
a few 10 32 ergs -1 . Both X-ray temperatures and luminosities are 
at variance with those displayed by normal Be stars. The origin of 
the X-ray emission in these objects is debated, it can either be as¬ 
sociated with accretion onto a white dwarf companion or due to 
magnetic int eraction between the star and the inner part of the de¬ 
cretion disc l lSmith & Robinsonl lT999). However, a recent analysis 
of the long term optical and X-ray light curves of y-Cas strongly 
favour magnetic activity as the source of the hard X-ray emission 
dMotch et alj|2015t) . The X-ray luminosity and hard thermal spec¬ 
trum of Source #5 indicate that the star is a likely new y-Cas ana¬ 
logue. 

We conclude that the X-ray emission of sources Source #2 and 
Source #3 is consistent with that expected from radiative shocks 
in the massive winds of single WR or Ofpe stars or from a wind¬ 
colliding binary. The supergiant Source #3 may also qualify as a 
SFXT in quiescence and the giant Be Source #5 is consistent with 
the class of y-Cas analogues. 

Source #1 and Source #4 were classified as M2I stars based 
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on the IR spectra taken at the WHT of 2MASSJ 17434750- 
2923098 and 2MASSJ 18521006+0012073 respectively. Red su¬ 
pergiant stars are descendents of OB stars that went through epis¬ 
odes of strong mass loss during their evolution on the giant branch. 
Although so far no RSG has been confirmed in X-rays, SctX-1 
has been su gges ted as a possib le HMXBs with a red supergiant 
companion iKaplan et al]l2007l) . This system has been classified 
as a likely symbiotic binary. Symbiotic binaries are wide binaries 
in which a white dwarf or a neutron star accretes matter from the 
wind of a companion red giant star (note that SctX-1 may have 
either a late giant or a supergiant luminosity class). They have 
soft X-ray spectra, and show optical emission line spectra. Re¬ 
cent studies, based on Swift/XRT observations, have nevertheless 
shown that this c lass of objects can also have hard X-ray spectra 
flLuna et alj|2013l) . Although we did not detect any emission line 
in the infrared spectra of Source #1 and Source #4 which would 
rule out a possible symbiotic nature of th e source, there are a few 
symbiotic binaries lacking emission lines dMunari & Z witte J 20021 : 
iHvnes et al.ll2013l) . The nature of the X-ray emission of Source #1 
and Source #4 remains unclear. Although they could belong to the 
class of HMXB with a late supergiant donor stars, it is more likely 
that 2MASSJ 17434750-2923098 and2MASS 18521006+0012073 
are spurious infrared identifications of Source #1 and Source #4. 

6.2 The nature of the hard X-ray sources 

We found 690 3XMM-2MASS-GLIMPSE unique associations 
with reliable X-ray parameters which allowed us to discern 
between soft and hard X-ray sources. Making use of SIMBAD 
and VizieR we compiled the list of previously classified objects 
(see Section 1231 The fraction of identified sources in the soft 
and hard X-ray bands are similar (18% versus 16%). In Table [6] 
we list the number of classified X-ray sources per class of object. 
Soft X-ray emitters are dominated by stars with spectral types O- 
M, representing about 90% of identified soft X-ray sources. Hard 
X-ray sources are dominated by HMXBs, 41% of the classified 
hard X-ray sources. Nevertheless the population of Wolf-Rayet 
and O-B stars being discovered in the hard X-ray band is con¬ 
stantly increasing. In this study we add three more sources to this 
new population (see Fig. rm About 81% of the identified hard X- 
ray sources are high-mass stars either single or in binary systems 
(WR+Be+HMXB+O/B from Table 0. 

However, the nature of six hard sources identified from the 
literature remains unclear. Source 3XMMJ182833.7-103702 was 
classified as a LMXB due to the occurence of type I X-ray bursts 
dComelisse et al J [2002t> . The X-ray burst showed low brightness, 
and furthermore it did not show a clear exponential cooling de¬ 
cay as can be seen from their figure 1. Since so far the spec¬ 
tral type of the secondary star is not confirmed we still con¬ 
sider this source as a LMXB candidate (note that this is not the 
soft LMXB presented in Table 0. Sources 3XMM J 155033.7- 
540916 and 3XMMJ183255.1-100741 have been classified as 
AGB candidates and 3XMM J184841.0-023929 as a YSO candid¬ 
ate bv lRobitaille et alj j 20081) . Sources 3XMM J182746.3-120605 
and 3XMM J182814.3-103728 have been classified as candidate 
HMXBs bv lMotch et aD J201oh (XGPS 10 and XGPS 15 respect¬ 
ively). The hard sample source 3XMMJ182746.3-120605 is the 
same as XGPS-10 in iMotch et al.l J2010I) . We note that the pos¬ 
ition used for XGPS-10 in the latter work (RA = 18:27:45.49 , 
Dec = -12:06:06.7) was derived from an early XMM-Newton re¬ 
duction pipeline in 2001 and is very significantly offset by 12.2 
arcsec with respect to that now provided by the 3XMM catalogue. 


Table 6. Statistics of source type. 



Total 

WR 

Be 

HMXB 

LMXB 

Nova O/B 

A-M 

AGN 

Unid 

Hard 

173 

5 

4 

11 

0 

1 

2 

3 

1 

146 



3% 

2% 

6% 


1% 

1% 

2% 

1% 

84% 

Soft 

517 

7 

0 

0 

1 

0 

30 

55 

0 

424 



1% 

0% 

0% 

0 

0% 

6% 

11% 

0% 

82% 

All 

690 

12 

4 

11 

1 

1 

32 

58 

1 

570 



2% 

<1% 

2% 

<1% 

<1% 

5% 

8% 

<1% 

83% 


Notes. The table shows the fraction of each type of source among a sample 
of 640 X-ray sources with infrared counterparts. The cross-match procedure 
and the selection criteria applied to this sources are presented in Sections [2.2l 
and 12.31 


Accordingly, the faint possible candidate deemed as uncertain by 
IMotch et all | |2010|) cannot be the actual counterpart. Our work 
shows that XGPS-10 is in fact identified with the nearby bright 
K=11.94 infrared source. To conclude, since the nature of these six 
sources is still not confirmed we prefer to leave them as unidenti¬ 
fied (i. e., among the 146 Hard Unid in Table 0, meaning that a 
confirmation of their nature is needed. 

Finally, three stars classified as active coronae have hard X-ray 
spectra. These three stars have all G spectral type. Although active 
stellar coronae have typically soft X-ray spectra, many BY Dra- 
conis and RS CVn binaries, i. e. main sequence and evolved bin¬ 
aries, have been detected in hard X-rays (e. g. UPeg lOsten et al J 
120071) . The three G stars in this sample could belong to these types 
of objects. 

We show the K s magnitude versus X-ray-to-infrared flux ra¬ 
tio log(f x /fK s ) for hard sources in Fig. |TT] highlighting the posi¬ 
tion of the different types of objects found (HMXBs, WR, Be stars 
and late-type stars). In our hard X-ray source sample, Wolf-Rayet 
stars dominate at X-ray-to-infrared flux ratios log(f x /f Ks ) < -2.6, 
Be stars cover a very narrow range of values (although there are 
only a few of these systems), with mean log(f x /f Ks ) ~ -2.6, and 
HMXB dominate at higher values, log(f x /f Ks ) > -2.7. It should 
be stressed that the sources we selected for identification at the tele¬ 
scope in our pilot study are in no way representative of the entire set 
of pre-selected hard X-ray sources. First, we selected the brightest 
infrared sources for reasons of efficiency. Since our X-ray sample 
is flux limited this implied selecting sources at the low end of the 
X-ray-to-infra-red flux ratio distribution. In addition, as shown in 
Fig. 0 our pilot sample lies among the most reddened or intrins¬ 
ically red stars. This is due in part to the presence of two possible 
M2I interlopers and to the two Wolf-Rayet stars. Finally, due to the 
location of the William Herschel Telescope, we could only survey 
the Northern part of the Galactic plane covered by GLIMPSE. 

More information on the likely source content of the hard 
sample can be derived from identifications made through archival 
catalogues and from the literature. Table[7]hsts HMXBs and y-Cas 
analogues present in our sample. 

We found that the various classes of sources were best separ¬ 
ated in a 3-D space whose axis are an infrared colour (giving in¬ 
formation on both intrinsic colours and amount of reddening) HR3 
(which is the best defined proxy of the X-ray spectral hardness) and 
the X-ray-to-infrared flux ratio (which provides information on the 
emission mechanism). The top panel in Fig.ll2lshows the distribu¬ 
tion of the hard sample in this 3-D space, while introducing an IR 
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Table 7. HMXBs and Be stars identified from the literature. 


Source Name 

Sp Type 

log(f x /f Ks ) 

HR3 

Ks 

Note 

HD 119682 

B0.5Ve 

-2.56 

-0.21 

6.95 

y-Cas analogue 

SS 397 

B0.5Ve 

-2.60 

+0.08 

8.27 

y-Cas analogue 

NGC 6649 WL 9 

Bl-1.5IIIe 

-2.73 

+0.15 

7.79 

y-Cas analogue 

IGR J16318-4848 

SgB[e] 

-2.60 

+0.89 

7.19 

supergiant X-ray binary 

EXO 1722-363 

BO-lIa 

-0.53 

+0.96 

10.67 

wind fed Sg accreting neutron star 

Set X-l 

late Sg 

-2.67 

+0.92 

6.55 

supergiant X-ray binary 

AX J1910.7+0917 

B 

+0.04 

+0.82 

11.81 

early B-type supergiant star 

AX 11749.2-2725+ 

B3 

-2.04 

0.976 

10.07 

wind-fed accreting neutron stars 

IGR J16465-4507 

09.51a 

-1.08 

+0.51 

9.84 

supergiant fast X-ray transient 

XTEJ1739-302 

08.5Iab 

-2.55 

+0.62 

7.43 

supergiant fast X-ray transient 

IGR J17544-2619 

091b 

-2.26 

+0.50 

8.02 

supergiant fast X-ray transient 

IGR J18410-0535 

Bllb 

-1.26 

+0.93 

8.93 

supergiant fast X-ray transient 

IGR 118450-0435 

091a 

-0.82 

+0.59 

8.95 

supergiant fast X-ray transient 

IGR 118483-0311 

B0.5Ia 

-2.59 

+0.82 

8.46 

supergiant fast X-ray transient 


Notes. ' The infrared counterpart of AX J1749.2-2725 from whi ch a spectral type B was determined has infrared magnitudes fainter than those we determine 
in this paper, implying the infrared counterpart is different. iKaur et all <2010 1 derived the X-ray position from XMM and Chandra observations and infrared 
imaging and spectroscopy was carried out at the NTT, e.g. with a better resolution than 2MASS. We note that our identification of the true counterpart of the 
source might be wrong. 


magnitude allow us to sort out identified sources from the rest of 
the unidentified population (see the bottom panel in Fig. 112b . 

The group of four y-Cas analogues (three previously iden¬ 
tified) share similar X-ray and infrared colours and exhibit very 
comparable average X-ray-to-infra-red flux ratios. As illustrated in 
the two panels in Fig. ED many hard sources have properties sim¬ 
ilar to those of y-Cas analogues, albeit with fainter magnitudes and 
expected larger distances and may indeed be good candidates for 
belonging to this particular class. 

Very few active coronae have hard enough X-ray spectra to 
qualify as a hard X-ray source. However, several of the hard sources 
with HR3 < 0.6, log(f x /fii s ) in the range of -2 to -3 and no infra-red 
excess could still be identified with extreme active coronae. 

Compared to y-Cas analogues, Wolf-Rayet stars display on 
average softer X-ray spectra. Consequently, several of the WR stars 
identified in the literature actually failed to qualify as hard sources 
due to their too soft HR4, although their HR3 were fulfilling our 
conditions. It is thus unlikely that many WR stars appear in our 
selection of hard sources with HR3 < 0.2. However, four WR stars, 
two identified in the literature and two from our work, are found 
among our sample of hard sources with HR3 > 0.2. Accordingly, 
we expect to discover new X-ray selected WR stars in the region 
comprising the hardest X-ray sources with low to medium X-ray- 
to-infrared flux rat ios. We note that the WRs in our sample fall 
in lMauerhan et all d20ld) best selection areas in the (K s -8.0/rm) 
colour (see their Figure 8). 

The top panel in Fig. [T2] shows the presence of several hard 
sources in the region bounded by HR3 > 0.4, log(f x /fK s ) > -3 and 
low to high infrared excess. This area is populated by supergiant 
X-ray binaries and prospects are therefore relatively good that we 
find new massive X-ray binaries in this sample. Their relatively 
modest X-ray-to-infrared flux ratios would rather suggest SFXTs 
or classical Be/X-ray transients in the close to quiescence state. We 
note that among our relatively small sample of high-mass X-ray 
binaries SFXTs seem to display bluer IR colours. 

The nature of the population of faint (Ks >11) sources with 
log(f x /f Xs ) in the range of -2 to -1 and HR3 < 0.5 (yellow spots 


in Fig. remains essentially unconstrained since only one iden¬ 
tified source shares these properties, an AGN. Our identified AGN 
has log(f x /f Ks ) = -1.90, HR3 = 0.60, Ks = 12.65 and H-4.5yu = 
1.88. We determined the t otal galactic extinction in the line of sight 
using the Schlegel maps dSchleeel et al.lfl998l) . and found no rela¬ 
tion between the measured HR2 or HR3 and the total N H for hard 
sources with K s > 11 and log(f x /fK s ) > -2.5. The nature of 
this group of faint sources is thus not dominated by a population 
of background AGN. Likewise, the significantly higher X-ray-to- 
infrared flux ratio displayed by these faint sources compared to 
brighter ones suggests that they do not represent a reddened and 
remote population of sources identified at shorter distances. 


7 CONCLUSIONS 

We presented a study which aims at finding hard X-ray high mass 
stars with low- to intermediate- X-ray -luminosity. For that purpose 
we cross-correlated the 3XMM-DR4 catalogue with the 2MASS 
and GLIMPSE catalogues. We pre-selected 173 hard X-ray sources 
with bright infrared counterparts and with harder X-ray spectra 
(hardness ratio) than expected for a normal stellar corona. Among 
these 173 hard X-ray sources, 25 were previously classified: 3 
WRs, 3 Be stars, 11 HMXBs, 1 LMXB, 1 Nova, 5 stars and 1 
AGN. In this pilot study we carried out infrared spectroscopy of 
five new high-mass star candidates. We derived spectral types, and 
stellar parameters for these systems. We confirmed three sources 
as being high mass hard X-ray stars with spectral types WN7- 
8h, Ofpe/WN9 and Be and with X-ray luminosities in the range 
~ 10 32 - 10 33 erg s _1 (0.2 - 12 keV). One source is compatible with 
a colliding-wind binary, another source is compatible with both a 
colliding-wind binary and a supergiant fast X-ray transient and one 
source is a y-Cas analogue candidate. We classified two 2MASS 
sources as supergiant stars, both with a M21 spectral type, and 
concluded that they are likely not the true counterpart of the X- 
ray sources, whose nature is still unkown, although they could also 
belong to the class of HMXBs with a late supergiant star. Finally, 
the distribution of hard X-ray sources in the parameter space made 
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Figure 12. Colours of hard X-ray sources. In both figures symbols are colour coded according to the ratio of X-ray-to-infrared fluxes. Unclassified sources 
are shown as filled circles, WRs as open circles (we also included soft WRs in these two plots), HMXBs as triangles. Be stars as filled squares and other stars 
with open squares. The position of the only AGN is shown with an upside down triangle. Targets observed in this study are highlighted with a five pointed star 
symbol. 
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WR Be HMXB LMXB Nova Stars AGN 


Figure 13. Fraction of soft sources (light grey) and hard sources (grey) per 
object type. In black we show the position of the three X-ray sources for 
which we give a classification in this study. 

of X-ray hardness ratio, infrared colours and X-ray-to-infrared flux 
ratio suggests that many of the unidentified sources are new y-Cas 
analogues, WRs and low Lx HMXBs. However, the nature of the 
X-ray population appearing at K s > 11 and exhibiting average X- 
ray-to-infrared flux ratios remains unconstrained. 


ACKNOWLEDGEMENTS 

We would like to thank our referee I. Negueruela for his use¬ 
ful comments which significantly helped improving the quality of 
this work. We also thank P. A. Crowther and C. K. Rosslowe for 
sharing with us their new calibrated infrared Wolf-Rayet absolute 
magnitudes. We acknowledge financial support from the ARCHES 
project (7th Framework of the European Union, n 313146). FJC 
acknowledges financial support from the Spanish Ministerio de 
Econofa y Competitividad under project AYA2012-31447. This 
publication makes use of data products front the Two Micron All 
Sky Survey, which is a joint project of the University of Massachu¬ 
setts and the Infrared Processing and Analysis Center/California 
Institute of Technology, funded by the National Aeronautics and 
Space Administration and the National Science Foundation. This 
work is based in part on observations made with the Spitzer Space 
Telescope, which is operated by the Jet Propulsion Laboratory. 
California Institute of Technology under a contract with NASA. 
This research has made use of NASA’s Astrophysics Data System. 
This research has made use of the SIMBAD and VizieR databases, 
operated at CDS, Strasbourg, France. Based on spectral-fitting res¬ 
ults from the XMM-Newton spectral-fit database, an ESA PRO- 
DEX funded project, based in turn on observations obtained with 
XMM-Newton, an ESA science mission with instruments and con¬ 
tributions directly funded by ESA Member States and NASA. 


References 

Anderson G. E. et al.. 2011, ApJ, 727, 105 
Anderson L. D., Bania T. M., Jackson J. M., Clemens D. P.. Heyer 
M„ Simon R„ Shah R. Y„ Rathborne J. M„ 2009, ApJS, 181, 255 
Becker R. H., White R. L., Helfand D. J., Zoonematkermani S., 
1994, ApJS, 91, 347 

Benjamin R. A. et al„ 2003, PASP, 115, 953 


Blum R. D„ Damineli A., Conti P. S„ 2001, AJ, 121, 3149 
Bonnarel F. et al„ 2000, A&AS. 143, 33 
Cardelli J. A., Clayton G. C., Mathis J. S., 1989, ApJ, 345, 245 
Carey S. J. et ah, 2009, PASP, 121, 76 
Castelli F., Kurucz R. L., 2004, ArXiv Astrophysics e-prints 
Chaty S., 2011, in Astronomical Society of the Pacific Confer¬ 
ence Series, Vol. 447, Evolution of Compact Binaries, Schmidto- 
breick L., Schreiber M. R., Tappert C., eds., p. 29 
Churchwell E. et ah, 2009, PASP, 121, 213 
Churchwell E. et ah, 2006, ApJ, 649, 759 
Clark J. S., Steele I. A., 2000, A&AS, 141, 65 
Cornelisse R. et ah, 2002, A&A, 392, 885 
Corral A., Georgantopoulos I., Rosen S., Watson M. G., Page K., 
Stewart G. C., 2014, ArXiv e-prints 
Covey K. R. et ah, 2007, AJ, 134, 2398 

Crowther P. A., Hadfield L. J., Clark J. S., Negueruela I., Vacca 
W. D„ 2006, MNRAS, 372, 1407 
Cutri R. M. et ah, 2003, VizieR Online Data Catalog. 2246, 0 
Cutri R. M. e., 2012, VizieR Online Data Catalog, 2311, 0 
Deharveng L. et al., 2010, A&A, 523, A6 
Deharveng L. et al., 2012, A&A, 546, A74 
Dempsey R. C., Linsky J. L., Schmitt J. H. M. M.. Fleming T. A., 
1993, ApJ, 413, 333 

den Herder J. W. et ah, 2001, A&A, 365, L7 

Ebisawa K. et ah, 2005, ApJ, 635, 214 

Egan M. P. et al., 2003. VizieR Online Data Catalog, 5114, 0 

Figer D. F., McLean I. S., Najarro F., 1997, ApJ, 486, 420 

Forster Schreiber N. M., 2000, AJ, 120, 2089 

Gagne M. et ah, 2011, ApJS, 194, 5 

Gagne M., Oksala M. E., Cohen D. H., Tonnesen S. K.. ud-Doula 
A., Owocki S. P, Townsend R. H. D., MacFarlane J. J., 2005, 
ApJ, 628, 986 

Gilfanov M„ 2004, MNRAS, 349, 146 
Grimm H.-J., Gilfanov M.. Sunyaev R„ 2002, A&A, 391, 923 
Giidel M„ Guinan E. F„ Skinner S. L„ 1997, ApJ, 483, 947 
Hamann W.-R„ Grafener G., 2004, A&A, 427, 697 
Hands A. D. P„ Warwick R. S., Watson M. G„ Helfand D. J„ 
2004, MNRAS, 351, 31 

Hanson M. M., Kudritzki R.-P, Kenworthy M. A., Puls J., Tokun- 
aga A. T„ 2005, ApJS, 161, 154 

Hanson M. M„ Rieke G. H„ Luhrnan K. L„ 1998, AJ, 116, 1915 
Helfand D. J„ Becker R. H„ White R. L„ Fallon A., Tuttle S„ 
2006, AJ, 131, 2525 

Hoare M. G. et ah, 2012, PASP. 124, 939 

Hofner P, Delgado H., Whitney B., Churchwell E., Linz H., 2002, 
ApJ, 579, L95 

Hong J. S., van den Berg M., Grindlay J. E., Laycock S., 2009, 
ApJ, 706, 223 

Horne K„ 1986, PASP, 98, 609 
Hynes R. I. et ah, 2013, ArXiv e-prints 
Indebetouw R. et ah, 2005, ApJ, 619, 931 
Jansen F. et ah, 2001, A&A, 365, LI 

Kaplan D. L., Levine A. M., Chakrabarty D., Morgan E. H., Erb 
D. K., Gaensler B. M., Moon D.-S., Cameron P. B., 2007, ApJ, 
661, 437 

Kaur R.. Wijnands R., Paul B., Patruno A.. Degenaar N., 2010, 
MNRAS, 402. 2388 

Kobulnicky H. A., Babler B. L., Alexander M. J., Meade M. R., 
Whitney B. A„ Churchwell E. B„ 2013, ApJS, 207, 9 
Levesque E. M., Massey P, Olsen K. A. G., Plez B., Josselin E., 
Maeder A., Meynet G., 2005. ApJ, 628, 973 




18 A. Nebot Gomez-Moran el al. 


Lopes de Oliveira R., 2007, PhD thesis, Institute de Astronomia, 
Geoflsica e Ciencias Atmosfericas, Universidade de Sao Paulo, 
R. do Matao 1226, 05508-090 Sao Paulo, Brazil Observatoire 
Astronomique, UMR 7550 CNRS, Universite Louis Pasteur, 11 
rue de TUniversite, 67000 Strasbourg, France 
Lucy L. B„ White R. L„ 1980, ApJ, 241, 300 
Luna G. J. M„ Sokoloski J. L„ Mukai K„ Nelson T„ 2013, A&A, 
559, A6 

Marigo R, Girardi L., Bressan A., Groenewegen M. A. T., Silva 
L„ Granato G. L„ 2008, A&A, 482, 883 
Marshall D. J., Robin A. C., Reyle C., Schultheis M., Picaud S., 
2006, A&A, 453, 635 
Martins F„ Plez B„ 2006, A&A. 457. 637 
Mateos S., Saxton R. D., Read A. M., Sembay S., 2009, A&A, 
496, 879 

Mauerhan J. C., Muno M. P, Morris M. R., Bauer F. E., 
Nishiyama S., Nagata T., 2009, ApJ, 703, 30 
Mauerhan J. C., Muno M. P., Morris M. R., Stolovy S. R., Cotera 
A., 2010, ApJ, 710, 706 

Messineo M., Zhu Q., Ivanov V. D., Figer D. F., Davies B., Men- 
ten K. M„ Kudritzki R. P„ Chen C.-H. R„ 2014, A&A, 571, A43 
Meynet G„ Maeder A., 2000, A&A, 361, 101 
Monet D. G. et al„ 2003, AJ. 125, 984 

Morris P. W., Eenens P. R. J., Hanson M. M., Conti P. S., Blum 
R. D„ 1996, ApJ, 470, 597 

Motch C., Guillout P., Haberl F., Pakull M., Pietsch W., Reinsch 
K„ 1997, A&A, 318, 111 

Motch C., Herent O., Guillout R, 2003, Astronomische Na- 
chrichten, 324, 61 

Motch C., Lopes de Oliveira R., Negueruela I., Haberl F., Janot- 
Pacheco E., 2007, in Astronomical Society of the Pacific Confer¬ 
ence Series, Vol. 361, Active OB-Stars: Laboratories for Stellare 
and Circumstellar Physics, Okazaki A. T., Owocki S. R, Stefl S., 
eds., p. 117 

Motch C., Lopes de Oliveira R., Smith M. A., 2015, ArXiv e- 
prints 

Motch C. et al., 2010, A&A, 523, A92 
Munari U., Zwitter T., 2002, A&A, 383, 188 
Muno M. P. et al., 2009, ApJS, 181, 110 
Nebot Gomez-Moran A. et al., 2013, A&A, 553, A12 
Negueruela I., Gonzalez-Femandez C., Marco A., Clark J. S., 
2011, A&A, 528, A59 

Negueruela I., Gonzalez-Fernandez C., Marco A., Clark J. S., 
Martmez-Nunez S., 2010, A&A, 513, A74 
Negueruela I., Smith D. M., Reig R, Chaty S., Torrejon J. M., 
2006, in ESA Special Publication, Vol. 604, The X-ray Universe 
2005, Wilson A., ed„ p. 165 

Nesterov V. V., Kuzmin A. V., Ashimbaeva N. T., Volchkov A. A., 
Roser S„ Bastian U„ 1995, A&AS, 110, 367 
Ochsenbein F., Bauer P., Marcout J., 2000, A&AS, 143, 23 
Oskinova L. M., Hamann W.-R., Feldmeier A., Ignace R., Chu 
Y.-H., 2009, ApJ, 693, L44 

Osten R. A., Drake S., Tueller J., Cummings J.. Perri M., Moretti 
A., Covino S„ 2007, ApJ, 654, 1052 
Peretto N„ Fuller G. A., 2009, A&A, 505, 405 
Pfahl E., Rappaport S., Podsiadlowski R, 2002, ApJ, 571, L37 
Pineau F.-X., Motch C., Carrera F., Della Ceca R.. Derriere S., 
Michel L„ Schwope A., Watson M. G„ 2011, A&A, 527, A126 
Predehl R, Schmitt J. H. M. M„ 1995, A&A, 293, 889 
Raguzova N. V., 2001, A&A, 367, 848 

Revnivtsev M., Sazonov S., Churazov E., Forman W., Vikhlinin 
A., Sunyaev R., 2009, Nature, 458, 1142 


Robitaille T. P. et al., 2008, AJ, 136, 2413 
Romano R, Ducci L., Mangano V., Esposito R, Bozzo E., Vercel- 
lone S„ 2014, A&A, 568, A55 
Rosslowe C. K„ Crowther P. A., 2015, MNRAS, 447, 2322 
Sanna A. et al., 2014, ApJ, 781, 108 

Sazonov S., Revnivtsev M., Gilfanov M., Churazov E., Sunyaev 
R., 2006, A&A, 450, 117 

Schlegel D. J., Finkbeiner D. R, Davis M., 1998, ApJ, 500, 525 
Schultheis M. et al„ 1999, A&A, 349, L69 
Sguera V. et al., 2006, ApJ, 646, 452 
Shara M. M. et al., 2009, AJ, 138, 402 

Skinner S. L., Zhekov S. A., Giidel M., Schmutz W., Sokal K. R., 
2010, AJ, 139, 825 

Skrutskie M. F. et al., 2006, AJ, 131, 1163 
Smith M. A., Lopes de Oliveira R., Motch C., 2015, in Astro¬ 
nomical Society of the Pacific Conference Series, Vol. 0, Bright 
Emissaries, Sigut A., C. J., eds., p. 0 
Smith M. A., Robinson R. D„ 1999, ApJ, 517, 866 
Sokal K. R„ Skinner S. L„ Zhekov S. A., Giidel M„ Schmutz W„ 
2010, ApJ, 715, 1327 

Spitzer Science C., 2009, VizieR Online Data Catalog, 2293, 0 
Stead J. J., Hoare M. G., 2011, MNRAS, 418, 2219 
Steele I. A., Clark J. S., 2001, A&A, 371, 643 
Striider L. et al., 2001, A&A, 365, L18 

Sugizaki M., Mitsuda K., Kaneda H., Matsuzaki K., Yamauchi S., 
Koyama K„ 2001, ApJS, 134, 77 
Tauris T. M., van den Heuvel E. P. J., 2006, Formation and evol¬ 
ution of compact stellar X-ray sources, Lewin W. H. G., van der 
Klis M., eds., pp. 623-665 

Taylor M., 2011, TOPCAT: Tool for Operations on Catalogues 
And Tables. Astrophysics Source Code Library 
Torrejon J. M., Schulz N. S., Nowak M. A., Testa R, Rodes J. J., 
2013, ApJ, 765, 13 

Turner M. J. L. et al., 2001, A&A, 365, L27 

van Belle G. T. et al., 1999, AJ, 117, 521 

Watson M. G. et al., 2009, A&A, 493, 339 

White R. L„ Becker R. H„ Helfand D. J., 2005, AJ, 130, 586 

Wright E. L. et al., 2010, AJ, 140, 1868 

Yamaguchi H., Ueno M., Koyama K., Bamba A., Yamauchi S., 
2004, PASJ, 56, 1059 


APPENDIX A: NOTES ON B AND BE HARD X-RAY 
STARS 

We briefly discuss the nature of the hard B and Be stars found in 
this study via cross-correlations with SIMBAD and our spectro¬ 
scopic observations. 

3XMM J134632.5-625524 The optical counterpart of 
3XMM J 134632.5-625524, HD 119682, was found to be a 
B0.5Ve star at a distance of 1.3 kpc and with X-ray luminosity 
L x ~ lO^ ergs -1 , belonging to the class of y-Cas analogues (see 
ITorreion et alj|2013l and references therein). It is a blue straggler 
belonging to the open cluster NGC 5281 at a distance of 1.3 kpc. 
3XMM J181046.7-194421 This X-ray source is likely associated 
to HD 312792 which was classified as B5 by INesterov et alj 
il995l) . Visual inspection of the optical spectra used for stellar 
classification was not possible. No estimates for the stellar lumin¬ 
osity class and distance to the star were found in the literature. 
3XMM J182746.3- 120605 aka XGPS- IJ182745-120606 is identi¬ 
fied with XGPS-10 jMotch et~aT]|201 (j) . We note that the position 
used for XGPS-10 in the latter work (RA = 18:27:45.49 , Dec = 
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-12:06:06.7) was derived from an early XMM-Newton reduction 
pipeline in 2001 and is very significantly offset by 12.2 arcsec 
with respect to that now provided by the 3XMM catalogue. 
Accordingly. the_faint possible candidate deemed as uncertain by 
iMotch et all l20ld ) cannot be the actual counterpart. Our work 
shows that XGPS-10 is identified with the nearby bright K= 11.94 
infrared source. 

3XMM J182814.3-103728 Based on the optical spectrum, which 
revealed Ha- emission and Paschen absorption lines, iMotch et all 
Eolch classified the source as a possible B star (source XGPS-15). 
Although the authors noted the possibility of the optical source 
being an interloper. At a distance of about lOkpc, consistent with 
the strong extinction observed for this source, its X-ray luminosity 
is ~ 2.4 x 10 34 ergs -1 . The source was classified by IMotch et all 
EoTch as a HMXB candidate, with a giant or supergiant star as 
secondary star. As for the case above, if the optical counterpart is 
a blue supergiant the source would not belong to the class of Be 
stars. 

3XMM J182828.0-102401 A spectrum taken for the optical 
count erpart USNO -A2.00750-13261865 revealed a B3 spectral 
type iMotch et alj 1201(3) lacking emission lines. No luminosity 
class and distance values are available in the literature. 

3XMM J183327.7-103524 Optical identification of the X-ray 
source 3XMMJ183327.7-103524 was firstly presented by 
IMotch etahl 12003 ). The optical spectrum of SS397 revealed 


a BOV star. Nebot Gomez-Moran et alj (120131) estimated 


distance to the source of 1.5 kpc and an X-ray luminosity of 
Lx = 4.4 x 10 32 ergs -1 , and suggested the source to be a likely 
y-Cas analogue. 

3XMM J183328.3-102408 Based on spectroscopic obser¬ 
vation^ of the optical counterpart USNO0750-13549725 C, 
iNebot Gomez-Moran et all 12013l) classified the source as a 
B1 -1.5 Ille star. They estimated a distance to the star of 1.7 kpc, 
a value consistent with that to the cluster NGC6649. USN00750- 
13549725 C is the brightest star in the cluster and it is a blue 
straggler. It’s X-ray luminosity is L x = 6 x 10 32 ergs -1 . The 
source belongs to the y-Cas analogue class. 

3XMM J190144.5+045914 We classified the source in this study 
as a Be star, at a distance in the range 1-1.8 kpc and a luminosity 

of L x =0.9 -3.0 x 10 32 erg s -1 . This source belongs to the class 

of y-Cas analogues. 

Out of the eight hard B stars presented in this section, two are 
Be/X-ray high mass binary candidates with X-ray luminosities 
L x ~ 10 34 erg s -1 (both with a giant or supergiant secondary star), 
four are y-Cas analogues with L x ~ 10 32 ergs -1 , and two are B 
stars with no emission lines detected so far. 





























